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Abstract 

Two  robust  model-based  controllers  and  two  decentralized  adaptive  controllers  are 
experimentally  evaluated.  Algorithm  evaluation  is  motivated  by  the  need  for  controllers 
with  good  high  speed  tracking  under  varying  payload  conditions.  The  test  case  is  a  PUMA- 
560  robotic  manipulator  operating  over  the  standard  test  suite. 

The  model-based  controllers  are  made  robust  by  the  addition  of  either  an  auxil¬ 
iary  input  term  or  an  adaptive  feedforward  compensator  based  on  Lyapunov  Theory.  The 
model-based  auxiliary  input  controller  (MBAIC)  adapts  the  gain  matrices  used  in  comput¬ 
ing  an  additional  torque  to  be  combined  with  model-based  feedforward  and  PD  feedback 
torques.  The  adaptive  model-based  controllers  adapt  the  assessment  of  the  manipulator 
parameters  used  in  calculating  feedforward  torque. 

The  decentralized  adaptive  controllers  are  based  on  either  Lyapunov  stability  or 
Popov  hyperstability.  These  controllers  calculate  feedforward,  feedback,  and  auxiliary 
torques  based  on  trajectory  errors  and  desired  trajectory  parameters.  The  gain  matrices 
used  to  multiply  these  quantities  are  adapted.  These  auxiliary  torque  components  are 
identical  to  those  used  in  the  MBAICs. 

Experimental  evaluation  provides  insight  into  the  potential  and  limitations  of  each 
method.  The  decentralized  digital  adaptive  control  algorithms  produce  an  unsatisfactory 
tracking  response.  Both  model-based  control  techniques  improve  the  manipulator’s  track¬ 
ing  response. 
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INVESTIGATION  OF  ADAPTIVE  CONTROLLERS 
FOR  PUMA  TRAJECTORY  TRACKING 

i 

I 


L  Introduction 


LI  Motivation 

The  United  States  Air  Force  proposes  using  robotic  manipulators  to  replace  humans  in 
hazardous  or  hostile  environments.  Robotic  telepresence  and  aircraft  refueling  applications 
will  require  a  manipulator  with  the  ability  to  emulate  the  performance  of^a  human  arm. 
These  manipulators  must  be  able  to  track  high-speed  trajectories  under  a  vatiety  of  payload 
conditions.  Research  on  alternative  robotic  control  methods  must  be  expdnded.  Current 
control  methods  need  to  be  more  fully  explored  and  tuned  to  derive  optimum  performance 
from  specific  manipulators.  1 

1.2  Problem  Statement 

A  robot  is  commanded  to  follow  a  desired  trajectory  by  its  control  algorithm.  Current 

I 

industrial  applications  use  controllers  employing  Proportional  Derivative  (PD)  or  Propor¬ 
tional  Integral  Derivative  (PID)  methods  for  applications  where  payloads  and  trajectories 
are  fixed  [6].  These  controllers,  which  feed  back  the  measured  errors  in  an  effort  to  com¬ 
pensate  for  disturbances  produced  by  unmodeled  robot  dynamics,  are  designed  based  on 
the  assumption  that  the  robot  dynamics  are  modeled  by  a  series  of  linear  second-order 
systems  [8],  PD  and  PID  controllers  perform  adequately  when  the  distuibances  are  small. 
Robot’s  dynamics  models  can  not  be  assumed  linear  under  conditions  of  variable  manip¬ 
ulator  speed.  Varying  payloads  violate  the  assumption  of  a  constant  model  [15].  Under 
either  of  these  two  conditions,  precomputed  constants  for  PD  or  PID  controllers  become 
invalid  and  cannot  provide  the  desired  level  of  trajectory  tracking  accuracy  [15]. 
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1.3  Research  Objective 

An  area  of  active  research  is  the  development  of  adaptive  control  algorithms  that 
account  for  the  nonlinear  dynamics  inherent  in  a  robotic  manipulator  [12,  13,  17,  19].  A 
basic  form  of  adaptive  control  is  the  model-based  approach.  In  model-based  control,  the 
controller  adapts  to  changes  in  robot  configurations  (speed  and  trajectory).  Experimen¬ 
tal  evaluations  have  demonstrated  the  potential  for  model-based  techniques  to  improve 
high-speed  trajectory  tracking  accuracy  [1,  7,  9].  Unknown  payload  variations  reduce  the 
effects  of  these  control  schemes  [9].  Control  algorithms  that  incorporate  knowledge  of  ma¬ 
nipulator  system  dynamics  and  are  robust  and/or  adapt  to  variations  in  those  dynamics 
caused  by  model  inaccuracies,  payload  variation,  and  environmental  interaction  will  be 
required  [11].  Control  algorithms  that  are  independent  of  the  dynamics  model  are  also  un¬ 
der  consideration  [20,  22,  23].  These  specific  algorithms  are  decentralized  or  independent 
of  joint-to-joint  interaction.  The  objective  of  this  research  is  to  continue  evaluation  of  an 
adaptive  model-based  control  algorithm  and  to  investigate  decentralized  digital  adaptive 
control  techniques. 

1-4  Method  of  Approach 

As  part  of  the  Air  Force  Institute  of  Technology  (AFIT)  Robotics  Laboratory’s  Gross 
Motion  Control  Program,  this  thesis  investigated  three  different  types  of  control  algorithms: 
Adaptive  Model-Based  Control  (AMBC),  decentralized  adaptive  control,  and  Model-Based 
Auxiliary  Input  Control  (MBAIC). 

The  comparisons  and  experimental  evaluations  performed  during  this  thesis  were 
conducted  under  the  AFIT  Robotic  Control  A'gorithm  Development  and  Evaluation  (AR¬ 
CADE)  environment..  ARCADE  is  hosted  on  a  VAXstation  III  and  has  both  serial  and 
parallel  connections  to  the  PUMA  computer  bus.  The  ARCADE  environment  is  discussed 
further  in  Chapter  3.  The  PUMA-560  manipulator  environment  is  chosen  because  a.ny  al¬ 
gorithm  that  performs  well  on  a  PUMA  should  work  even  better  on  a  more  modern  design 
[11].  Plots  of  the  trajectory  errors  will  be  produced  to  display  the  error  data  collected. 
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The  AMBC  controller  is  based  on  the  robot’s  dynamics  model.  The  strength  of 
an  AMBC  controller  is  its  ability  to  alter  the  robot’s  manipulator  parameters  to  improve 
tracking  response.  Repeating  the  same  trajectory  several  times,  and  allowing  the  controller 
to  retain  previously  computed  parameters,  is  referred  to  as  learning.  A  nominal  approx¬ 
imation  of  these  parameters  may  be  used  to  start  the  controller  (the  controller  would  be 
initialized),  or  the  parameters  could  be  set  to  zero  (uninitialized).  The  AMBC  controller 
had  already  been  implemented  and  evaluated  in-house  [12,  13].  This  software  was  modified 
to  include  static  friction  parameters  in  an  attempt  to  improve  initial  and  endpoint  errors. 
The  dynamics  model  was  rearranged  to  test  the  importance  of  friction  compensation.  All 
three  implementations  of  the  AMBC  controller  were  evaluated  for  single-run  trajectories 
and  learning  tests  using  initialized  and  uninitialized  parameters. 

Decentralized  adaptive  control  algorithms  are  independent  of  the  robot  dynamics 
model.  All  matrices  used  in  calculating  control  torques  are  diagonal  to  eliminate  inter¬ 
joint  dependence  (thus  the  term  “decentralized”).  Seraji’s  Lyapunov-based  decentralized 
digital  adaptive  controller  was  implemented  by  Leahy  [20,  10].  New  software  subroutines, 
based  loosely  on  this  code,  were  written  to  test  Tarokh’s  hyperstability-based  algorithm. 
These  controllers  use  the  products  of  adaptive  matrices  and  desired  trajectory  components 
or  position  errors  to  calculate  feedback,  feedforward,  and  auxiliary  torques.  After  tuning, 
Tarokh’s  controller  was  evaluated  for  trajectory  and  payload  independence.  The  test  results 
were  also  compared  to  Leahy’s  work  on  Seraji’s  controller  [10]. 

Model-Based  Auxiliary  Input  Controllers  combine  the  dynamics  model  of  the  AMBC 
algorithm  with  the  Auxiliary  Torque  of  decentralized  adaptive  controller.  The  model- 
based  feedforward  torque  is  computed  from  desired  trajectory  components.  The  controller 
uses  PD  feedback  torques.  Leahy  implemented  and  tested  a  Seraji-Based  MBAIC  [10].  A 
Tarokh-based  MBAIC  was  constructed,  tuned,  tested,  and  compared  to  Leahy’s  work. 

All  tests  were  conducted  over  a  proposed  set  of  standard  trajectories  that  excite  all 
the  PUMA  dynamics  and  allow  a  rigorous  test  of  algorithm  performance  [13].  All  of  the 
trajectories  were  traversed  under  no-load  conditions.  Motor  saturation  constraints  limited 
the  range  of  payload  testing  available.  Some  of  the  test  trajectories  were  run  with  a  2 
kilogram  brass  disk  attached  to  the  robot’s  end  effector.  The  mounting  of  this  simple  pay- 
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load  on  the  end-effector  was  an  adequate  approximation  of  the  very  real  task  of  robot  tool 
changing  and  covered  a  large  subset  of  the  payload  adaptation  tasks  required  for  a  robot  to 
emulate  human  arm  performance.  The  payload  was  heavy  enough  to  produce  a  significant 
degradation  in  algorithm  performance  without  driving  the  motors  into  saturation  [10]. 

1.5  Accomplishments 

Evaluations  of  an  existing  16-parameter  Adaptive  Model-Based  Controller  (AMBC) 
[13]  demonstrated  that  the  addition  of  three  static  friction  parameters  would  not  improve 
system  response  to  a  wide  variety  of  trajectory  and  payload  conditions.  A  13-parameter 
AMBC  controller  that  concentrated  heavily  on  friction  compensation  did  not  provide  suf¬ 
ficient  tracking  capability  but  showed  the  inability  of  an  AMBC  to  maintain  tracking 
accuracy  when  it  adapts  parameters  with  incorrect  signs. 

Examination  of  two  decentralized  adaptive  control  algorithms  [20,  22]  revealed  their 
limitations  and  unsuitability  for  gross  motion  control  applications.  This  was  AFIT’s  initial 
examination  of  Tarokh’s  decentralized  adaptive  controller  and  the  second  glance  at  Seraji’s. 
However,  Model-Based  Auxiliary  Input  Controllers  showed  promise.  These  controllers  are 
based  on  the  auxiliary  torque  components  derived  from  the  decentralized  controllers  and 
an  SMBC  controller.  The  MBAIC  controller  is  simple  to  implement,  relatively  easy  to 
tune,  and  provide  tracking  comparable  to  AMBC  controllers.  The  MBAIC  needs  to  be 
examined  further  to  investigate  potential  uses. 

1.6  Organization 

The  remainder  of  this  thesis  is  presented  in  five  chapters.  Chapter  2  presents  a  litera¬ 
ture  review  of  adaptive  model-based  control  schemes  and  decentralized  adaptive  control  al¬ 
gorithms  currently  under  investigation  at  the  Air  Force  Institute  of  Technology.  Chapter  3 
discusses  the  tuning,  testing,  and  analysis  of  the  adaptive  model-based  control  algorithm. 
Chaptei  4  coveis  the  tuning,  testing,  and  analysis  of  the  decentialized  adaptive  control  al¬ 
gorithms.  Chapter  5  covers  the  tuning,  testing,  and  analysis  of  the  Model-Based  Auxiliary 
Input  Control  algorithms  Finally,  Chapter  6  presents  conclusions  and  recommendations 
for  further  research. 


II.  Literature  Review 


2.1  Introduction 

This  chapter  is  divided  into  two  main  sections.  First,  model-based  robotic  control  is 
discussed  as  a  preiude  to  examining  an  adaptive  model-based  control  algorithm  that  has 
demonstrated  its  effectiveness  in  previous  in-house  evaluations  [11,  12,  13].  The  second 
portion  of  this  chapter  explores  two  decentralized  adaptive  control  algorithms  [20,  22,  23]. 
One  controller  is  based  on  Lyapunov  stability  criteria  whil°  the  other  is  formulated  using 
Popov  hyperstability  theory.  These  decentralized  algorithms  are  independent  of  robot 
dynamics  models. 

2.2  Adaptive  Model-Based  Control  Algorithm 

2.2.1  Model-Based  Control.  A  dynamics  model  of  the  robot  manipulator  must  be 
developed  in  order  to  use  model-based  control  techniques.  Ignoring  motor  dynamics  the 
equations  of  motion  may  be  written  in  vector  form  as  [4,  18]: 

r(t)  =  H  (0)0  +  C{0, 0)6  +  G  (0)  (2.1) 


where: 

n  is  the  number  of  degrees  of  freedom  of  the  robot  manipulator 

6,0,0  are  n-dimensional  vectors  of  joint  positions,  velocities,  and  accelerations 

rm(t)  is  an  n-dimensional  vector  of  joint  motor  torques 

H(0)  is  the  n  X  n  symmetric  positive-definite  inertia  matrix 

C(0, 6)0  is  an  n-dimensional  vector  of  centripetal  and  Coriolis  torques 

G(0)  is  an  ?j-dimensional  vector  of  gravity  torques 

The  end  result  of  model-based  control  is  to  obtain  an  equation  of  the  form: 


r(i)  =  rff+Tfb 


(2.2) 
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with: 


Tjf  =  feedforward  compensation  torques 
r/i,  =  feedback  compensation  torques 

The  right-hand  side  of  Equation  (2.1)  develops  the  feedforward  torques  which  may 
be  expressed  by: 

r//  =  H(0)0  +  C(0,0)0  +  G(0)  (2.3) 

where  *  denotes  modeled  values.  The  feedforward  torque  may  also  be  denoted  by  a  re¬ 
gressor  formulation  which  defines  the  torque  as  a  function  of  a  linear  parameter  vector 
a: 

rfJ  =  Y(9,9,9)k  (2.4) 

The  H,  C,  and  G  terms  are  expressed  in  terms  of  the  parameter  vector  and  combined  to 
form  the  regressor  matrix  Y.  The  modeling  attributed  to  the  H,  C,  and  G  matrices  is 
transferred  to  the  manipulator  parameters. 

The  feedback  torques  are: 

Tjt  =  K„e  +  Kpe  (2.5) 

where 


e  is  the  joint  velocity  error,  9^-9 
e  is  the  joint  position  error,  9^-9 

9d,9d  are  n-dimensional  vectors  of  desired  joint  positions  and  velocities 
K„,KP  are  diagonal  nxn  matrices  of  velocity  and  position  feedback  loop  gains 

The  ideal  situation,  arrived  at  by  substituting  Equations  (2.3)  and  (2.5)  into  Equation 
(2.2),  equating  the  result  to  Equation  (2.1)  and  assuming  that  the  modeled  values  are  equal 
to  the  actual  values,  is  denoted  by: 


e  T  K^.e  T  =  0 


(2.6) 
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Equation  (2.6)  states  that,  with  the  appropriate  Kp  and  K„  values,  the  model-based  con¬ 
troller  uniformly  rejects  errors  over  the  manipulator’s  entire  trajectory.  This  equivalence 
assumption  of  modeled  and  actual  values  ignores  uncertainties  that  arise  from  time  varia¬ 
tions  in  the  robot’s  transfer  function  produced  by  the  operating  environment,  component 
variation  due  to  wear,  and  manufacturing  irregularities  [8]. 

2.2.2  Adaptive  Model-Based  Control.  Uncertainties  in  the  dynamics  model,  Equa¬ 
tion  (2.3),  are  passed  on  to  the  model-based  controller  as  disturbances.  The  controller 
can  maintain  good  high-speed  tracking  performance  in  the  presence  of  small  disturbances. 
Unknown  payload  parameters  cause  degradation  of  the  controller’s  performance.  Some 
adaptive  control  algorithms  estimate  full  dynamics  feedforward  torque  compensation,  with 
unknown  manipulator  and  payload  parameters,  on  line.  These  algorithms  may  be  supple¬ 
mented  by  PD  feedback  torques.  Adaptive  control  schemes  based  on  Lyapunov  stability 
theory  have  been  developed  by  Sadegh  and  Horowitz  [19],  Slotine  and  Li  [21],  and  others 
[17,  24], 

For  this  discussion,  motor  dynamics  were  ignored  to  facilitate  the  presentation.  Dur¬ 
ing  the  actual  implementation  they  must  be  considered.  After  defining  a  control  law  of: 

r  =  H(0)0d  +  C(M)0rf  +  G(0)  +  Kpe  +  Kl;e  (2.7) 

=  Y(0,0A,0d)a  +  Kpe  +  K„4  (2.8) 

Slotine  and  Li  derive  an  adaptive  control  algorithm  from  a  Lyapunov  stability  analysis 
[21].  They  used  a  Lyapunov  function  candidate  of: 

V(t)  =  ^(eTH(<9)e  +  aTra  +  eTKpe)  (2.9) 


where: 

T  is  a  symmetric  positive-definite  matrix  (usually  diagonal) 
a  is  the  vector  of  manipulator  parameter  estimates 
a  =  a(/.)  -  a  is  the  parameter  estimation  error  vector 


Differentiating  Equation  (2.9)  and  using  the  property  of  skew  symmetry  described  in 
[21]  to  eliminate  the  |e  (H  -  2C)e  term  yields: 

V(t)  =  eT[r  -  H(0)0rf  —  C(0, 9)0  d  -  G(0)  —  Kpe]  +  aTTd  (2.10) 

=  eT[r-  Y(6>,0,^,0rf)a-Kpe]  +  aTrl  (2.11) 

After  substituting  Equation  (2.7)  for  r  and  defining: 

H(0)  =  H(0)-H(0) 

6(0,0)  =  6(0,0) -C(0,0) 

G(0)  =  G(0)-G(0) 

the  derivative  equation  (Equation  (2.10))  becomes: 

V(t)  =  -eTKue-eT[H(6>)0<i  +  C(0,0)0(i  +  G(0)]  +  aTri  (2.12) 

=  —eTK„e  —  eTY(0, 0, 0J,  0j)a  +  aTra  (2.13) 

Transposing  all  the  elements  in  the  second  term  of  the  right-hand  side  of  Equa¬ 
tion  (2.13)  and  combining  the  last  two  terms  leads  to: 

V(t)  =  —eTKve  +  aT[r4  -  YT(0, 0, 0d,  0d)e]  (2.14) 

For  Lyapunov  stability  to  hold  V(t)  must  be  zero  or  negative  [21].  Knowing  that  K„ 
is  defined  as  positive  definite,  then  setting: 

rk-  YT(o,d,ed,ed)k  =  o  (2.15) 

leaves: 

V(t)  =  -eTKve  <  0  (2.16) 

and  the  system  is  stable. 
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Rewriting  Equation  (2.15)  to  solve  for  the  unknown  vector  a  results  in: 


4  =  r-1YT(0,0)0d,0d)e  (2.17) 

Since  a  is  constant,  a  =  a,  and  we  have  an  equation  for  determining  an  accurate  model  for 
T//- 


2.2.3  Experimentation  Algorithm.  The  model-based  algorithm  for  the  output  torque 
used  by  Leahy  and  Whalen  [13],  adapted  from  Sadegh  and  Horowitz  [19],  contains  an  aux¬ 
iliary  input,  raX}  plus  feedforward  and  feedback  terms.  It  is: 

r  -  rff  +  TJb  +  T  ax  (2.18) 

The  feedback  portion  of  the  control  law  is: 

TJb  =  K„[(0d  -  9)  +  A(0d  -  e)]  (2.19) 

where  A  is  a  constant  matrix  with  eigenvalues  in  the  right-half  complex  plane.  Feedback 
torque  is  related  to  Equation  (2.16),  which  was  Slotine  and  Li’s  expression  for  Lyapunov 
stability.  A  is  selected  to  be  Kp/Kv  in  order  to  implement  PD  feedback  as  expressed  in 
Equation  (2.5). 

The  feedforward  component  of  this  equation  is  broken  down  into  two  parts  to  allow 
for  precalculation  and  off-line  storage  of  the  feedforward  torques  that  are  functions  of 
the  known  parameters,  an,  and  desired  joint  positions,  velocities,  and  accelerations.  The 
divided  equation  is: 


tss  =  Y(ed,ed,ed)& 

=  Yi(6d,6d,  0d)a.  ■\-Y-2(6d,dd,  6d)an  (2.20) 

The  selection  of  Equation  (2.7)  as  the  control  law  allows  the  cancellation  of  terms 
associated  with  the  known  manipulator  parameters,  leaving  only  the  unknown  parameters 
to  be  estimated  by  a.  The  regressor  matrix  does  not  depend  on  the  actual  manipulator 


acceleration,  but  on  the  velocity  and  acceleration  of  the  desired  trajectory  [17],  which 
eliminates  sensitivity  problems  that  may  be  induced  by  imprecise  measurements  if  on-line 
measurements  of  manipulator  acceleration  were  available  [21]. 

The  regressor  matrix  differs  from  Slotine  and  Li’s  version  because  it  is  only  depen¬ 
dent  on  desired  trajectory  components,  not  a  combination  of  actual  and  desired  tra¬ 
jectory  attributes.  The  feedforward  torque  component  based  on  known  components, 
Y2(Qd,Qd>8d)&n)  may  be  precomputed  and  stored  off-line.  The  portion  of  the  regres¬ 
sor  associated  with  the  unknown  parameters,  Y* (Od,0d,&d),  can  also  be  precomputed. 
This  approach  leads  to  a  reduction  in  the  number  of  on-line  calculations  which  permits 
the  use  of  an  increased  sampling  rate  [16].  The  main  advantage  of  using  desired  trajectory 
quantities  is  that  the  noisy  measurements  of  actual  velocity  and  acceleration  are  avoided. 
Other  advantages  afforded  by  adaptive  model-based  algorithms  are  increased  robustness 
and  the  capability  of  the  manipulator  to  “learn”  in  the  case  of  repetitive  tasks. 

The  parameter  adaptation  for  evaluating  the  unknown  manipulator  parameters  on 
line  was  accomplished  through  the  use  of  a  Desired  Compensation  Adaptation  Law  (DCAL). 
Sadegh  and  Horowitz’s  simulations  [19]  and  experimental  results  [16]  show  that  the  use  of 
a  DCAL  allowed  their  manipulator  to  maintain  stability  and  retain  parameter  convergence 
in  the  presence  of  noise  and  a  constant  adaptation  signal.  The  adaptation  law  is: 

a  =  J  T-1YT(Od,ed,Od)[(dd-0)  +  A(ed-O))  (2.21) 

This  adaptation  equation  is  nearly  the  same  as  Equation  (2.17),  except  for  the  addition  of 
the  weighted  position  error  term  A (Od  —  6).  This  regressor  matrix,  Y,  which  contains  only 
desired  trajectory  components,  can  be  precalculated  and  stored  off-line.  As  in  the  case  of 
the  precomputed  feedforward  torques,  using  a  regressor  matrix  based  on  desired  trajectory 
values  enhances  robustness  and  avoids  noise  corruption  of  the  error  and  adaptation  signals. 

The  auxiliary  torque  term  provides  compensation  for  additional  errors  introduced  by 
modifying  the  adaptive  controller  [19].  The  auxiliary  torque,  when  used,  is  evaluated  by: 

Tax  =  Vn\e\7[(0d  -0)  +  A (dd  -  6)}  (2.22) 
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This  function  was  set  to  zero  [13]  since  the  small  errors  anticipated  from  using  this  algo¬ 
rithm  make  this  torque  component  insignificant. 

2.2.4  Friction  Models .  Several  authors  [19, 11, 17]  note  the  importance  of  friction  in 
the  manipulator  dynamics  model.  The  different  types  of  friction  are  viscous,  Coulomb,  and 
static  friction  [8].  Viscous  friction,  rv(t))  is  characterized  as  a  linear  relationship  between 
applied  force  and  velocity.  Coulomb  (or  running)  friction,  rc(2),  has  a  constant  amplitude 
with  its  sign  dependent  on  the  sign  of  the  velocity.  Static  friction,  r5(Z),  is  the  resistance  to 
initial  motion  which  vanishes  with  the  onset  of  movement.  The  sign  of  the  static  friction 
term  is  such  that  the  force  opposes  motion.  These  forces  are  defined  respectively  by  [8]: 


t„(/)  =  B9(t) 

(2.23) 

rc(t )  =  Tcsgn(6 ) 

(2.24) 

rs(t)  =  ±r,\is0 

(2.25) 

The  static  friction  term  has  also  been  expressed  as  a  constant  times  an  exponential 
term  [2,  3].  Canudas  de  Wit  and  Seront  use  Tustin’s  model  to  express  the  overall  friction 
model  as  [5]: 

r(0)  =  [ao  +  aie-T°^  +  ot2\0\ }sgn(0)  (2.26) 

The  terms  ao,  a  1,  and  a 2  correspond  to  the  rc,  r5,  and  D  respectively  from  the  previ¬ 
ous  definitions.  The  exponential  term  allows  for  the  static  friction  term  to  be  a  rapidly 
decreasing  function  rather  than  a  spike.  The  term  r0  is  the  lubricant  coefficient  which  is 
between  20  and  60  for  most  materials.  These  authors  suggest  using  40  for  the  lubrication 
coefficient. 

Friction  has  a  significant  effect  on  the  performance  on  many  robotic  manipulators, 
especially  those  with  gear  reduction.  Friction  compensation  is  a  challenge  because  its 
functional  dependence  on  the  joint  variables  is  difficult  to  model  [17].  Leahy  and  Whalen 
state  that  their  adaptive  algorithm  benefits  from  its  ability  to  adapt  Coulomb  and  viscous 
friction  terms  [12].  This  adaptation  is  shown  to  reduce  the  tracking  errors  on  both  ends 
of  their  test  trajectories  where  velocity  is  low  and  friction  and  gravity  terms  dominate  the 
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dynamics.  Additional  friction  compensation  compensation  in  the  low  velocity  regions  of  the 
trajectory  may  further  reduce  tracking  errors.  The  friction  component  that  is  unmodeled 
by  both  Sadegh  and  Horowitz  [19],  and  Leahy  and  Whalen  [13],  is  static  friction. 

2.3  Decentralized  Digital  Adaptive  Control 

2.3.1  Introduction.  Conventional  industrial  controllers  are  based  on  independent 
joint  controllers  where  each  individual  joint  is  controlled  by  a  simple  position  servo-loop 
with  predefined  gains.  This  method  is  inadequate  for  applications  involving  precise  tracking 
of  fast  trajectories  under  uncertain  payload  conditions  [10,  20].  Control  systems  proposed 
to  handle  these  operating  conditions  include  model-based  and  performance-based  meth¬ 
ods.  ?erformance-based  methods  do  not  involve  the  use  of  the  complex  mathematical 
manipul  rtor  dynamics  model  in  the  foimulation  of  the  control  law. 

2.3.2  Lyapunov- Based  Decentralized  Adaptive  Control.  Seraji  proposed  a  simple 
Lyapunov-based  decentralized  adaptive  independent  joint  control  technique  [20].  This 
technique  does  not  require  any  knowledge  of  the  manipulator  dynamics  model,  model 
parameters,  or  payload  parameters.. 

The  coupled  manipulator  dynamics  model  needs  to  be  massaged  to  develop  a  decen¬ 
tralized  control  scheme.  Seraji’s  dynamics  model  is  [20]: 

T  =  H(0)0  +  C(0, 0)0  +  G(0)  +  S(0)  (2.27) 

where: 

T (t)  is  the  7i  X  1  vector  of  applied  joint  torques 
S(0)  is  the  ft  x  1  frictional  torque  vector 

This  model  can  be  broken  down  into  scalar  differential  equations  of  the  form: 

l\(t)  =  +  c,(0, 0)e,  +  g  ,(0)  +  Si(0)  (2.28) 

;= 1 
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where  i  denotes  the  ith  element  through  the  nth  joint  and  the  lower  case  letters  h,c,g,s  are 
the  individual  joint  variables  of  the  quantities  represented  by  the  same  upper  case  letters. 
This  model  is  equivalent  to: 


(2.29) 


where: 


di(9, 9, 9)  =  £  MW*)  +  a(0,  0)6i  +  gi(6)  +  3,(0)  (2.30) 

J  = 1 

Now  each  joint  is  modeled  by  an  input-output  dynamic  equation  with  the  joint  torque 
Tx  as  the  input,  the  joint  angle  0t*  as  the  output,  and  d{  treated  as  a  disturbance  torque. 
From  the  system  viewpoint,  the  term  dt  represents  the  coupling  between  the  ith  joint  and 
other  joints  and  nonlinear  dynamics.  The  joint  control  laws  are  restricted  to  be  decentral¬ 
ized,  meaning  that  each  controller  operates  on  only  its  own  joint  with  no  interchange  of 
information  between  joints.  In  other  words,  the  decentralized  manipulator  control  prob¬ 
lem  is  to  design  a  set  of  independent  joint  controllers  in  which  the  zth  controller  generates 
the  joint  torque  Tt(t)  by  acting  only  on  the  joint  angle  trajectory  &i(t)  and  the  desired 
trajectory  0^(0,  and  ensures  that  6t(t)  tracks  0^(0  asymptotically.  Seraji’s  proposed  de¬ 
centralized  controller  is  [20]: 

T «(0  =  /*(0  +  [&to(0c*(0  ^  +  [®o(0®*(0  +  5ii +  9*2(0 9di(t)}  (2.31) 

where,  for  each  individual  joint: 

9di(t)  is  the  desired  trajectory, 

9t (t)  is  the  joint  angle  trajectory,  and 
ct(/)  is  the  joint  tracking  error,  0^,(0  -  9t(t). 

This  control  law  has  three  distinct  components: 

1.  The  auxiliary  signal,  (/),  partially  compensates  for  dt(t)  and  improves  tracking 
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2.  The  term  [fc,oV)e«(*)  +  ^>i(0^»(0]  *s  the  adaptive  position- velocity  feedback  controller. 
The  gains,  fc,o(t)  and  kn(t),  are  adjustable. 

3.  The  term  [<fro(0^t(0  +  9«x(0^«(0  +  Qi2(t)^di(t)]  is  the  adaptive  position-velocity- 

accelerator  feedforward  controller.  Its  gains,  9,o(t),  9,i(t),and  are  also  ad¬ 

justable. 


Equation  (2.31)  can  be  rewritten  as: 

m  =  w) + + £«,  «*$«) 

i= o  i=° 


(2.32) 


where  the  superscripts  denote  derivatives  with  respect  to  time.  The  controller  adaptation 
laws  are  based  on  the  weighted  error  rt(t)  as  follows: 
weighted  error: 

ri(t)  =  wviCi(t )  +  wmet(t)  (2.33) 


auxiliary  signal: 

fi(t)  =  /,( 0)  +  Si  f  ri(t)dt  +  pin(t ) 
JO 

feedback  gains: 


r  ^ 

=  k{j( 0)  +  a{j  /  ri(t)e\3\t)dt.  +  j  =  0, 

Jo 


(2.34) 


(2.35) 


feedforward  gains: 

=  90(0)  +  7 ij  f  ri(t)9(J)(t)dt  +  3  =  0,1,2  (2.36) 

Jo 

where  {<$t, are  any  positive  scalar  integral  adaptation  gains,  Au}  are  zero 

or  any  positive  scalar  proportional  adaptation  gains,  and  {wpiywm}  are  positive  scalar 
weighting  factors  demonstrating  the  relative  magnitudes  of  the  position  and  velocity  errors 
used  in  forming  the  weighted  error. 

This  decentralized  controller,  consisting  of  several  independent  joint  controllers,  is 
purported  to  have  advantages  over  a  single  centralized  controller  foi  the  entire  manipu- 
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lator  [20].  The  decentralized  controller  is  more  computationally  efficient.  For  an  n-joint 
manipulator,  the  decentralized  controller  computes  6n  adaptive  gains  each  sample  period 
versus  5 n2  +  n  adaptive  gain  computations  for  a  centralized  controller.  System  performance 
could  be  further  enhanced  by  employing  n  simple  and  fast  microprocessors  for  parallel  pro¬ 
cessing  and  distributed  computing.  A  second  attribute  of  this  control  scheme  is  reliability 
and  fault  tolerance.  Erroneous  joint  position  readings  from  a  joint  encoder  would  affect 
the  entire  structure  of  a  centralized  controller  but  only  one  joint  of  the  decentralized  sys¬ 
tem,  A  limitation  of  the  decentralized  controller  is  the  possibility  of  tracking  errors  caused 
by  the  absence  of  complete  inter-joint  coupling  compensation.  The  decentralized  adaptive 
control  algorithm  accepts  this  deficiency  as  prrt  of  the  trade-off  in  attaining  computational 
simplicity,  fault  tolerance,  and  robustness. 

Seraji’s  experiments  were  conducted  on  a  PUMA-560  industrial  robot  over  a  single 
three  joint  trajectory  from  (0°,0°,0°)  to  (60°, -60°,  60°)  in  3  seconds  [20].  His  exper¬ 
imental  results  did  not  address  the  issues  of  trajectory  or  payload  independence  using 
the  decentralized  control  algorithm.  Seraji’s  design  parameters,  contained  in  Table  2.1, 
reduced  the  control  torque  to: 

T  (t)  =  f  (t)  +  ko  +  ki(£)e(2)  (2.37) 

which  uses  only  the  auxiliary  and  feedback  torque  components  available.  The  proportional 
pieces  of  the  f,  ko,  and  ki  adaptive  matrices,  pr,  /30re ,  and  /3\re  respectively,  were  also 


Table  2.1.  Seraji’s  resign  Parameters 


Auxiliary  terms: 

6i  =  [50, 50,50) 

Pi  =  [0,0,0] 

Feedback  terms: 

a,o  =  [100,100,100] 
a,i  =  [800,800,800) 

Pio  =  [0,0,0] 

Pix  =  [0,0,0] 

Feedforward  terms: 

7*o  =  [0,0,0] 

7.1  =  [0,0,0] 

7.2  =  [0,0,0] 

A;o  =  [0,0,0] 

An  =  [0,0,0] 

A.2  =  [o,  o,  uj 

Weighted  error  terms: 

wvl  =  [20,20,4] 
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omitted.  Leahy  [10]  further  evaluated  this  algorithm  using  Seraji’s  design  parameters  and 
initial  conditions  over  a  trajectory  of  (-50°,  -135°,  135°)  to  (45°,  -85°,  30°)  in  1.5  seconds 
to  excite  all  of  the  PUMA’s  dynamics.  Manipulator  performance  was  degraded.  Leahy 
concluded  that  this  control  scheme’s  performance  is  heavily  dependent  on  accurate  initial 
compensation  of  static  dynamic  forces  and  the  position  gain  adaptation  properties  of  the 
algorithm  are  minimal. 

Tarokh  developed  another  dec?  ntralized  control  scheme  based  on  Popov  hyperstability 
theory.  Hyperstability  theory  is  used  to  enhance  design  flexibility.  Tarokh’s  controller  is 
discussed  in  the  next  section. 

2.3.3  Discrete-Time  Adaptive  Control.  Tarokh  developed  a  discrete-time  adaptive 
control  algorithm  employing  Popov  hyperstability  [22,  23].  The  dynamics  model,  similar 
to  Equation  (2.27),  is  written  as: 

u(<)  =  B2  (O,b,m)0  +  B  i(6,0,m)6  +  Bo  (0,0, m)0  (2.38) 

where  u(t)  is  the  nx  1  torque  vector  and  B2(-)>  Bi(-),  and  Bo(-)  are  matrices  with  elements 
that  are  nonlinear  functions  of  payload  mass,  m(t),  and  manipulator  joint  position  and 
velocity. 

The  dynamics  mode'  can  be  further  reduced  to  Equation  (2.39)  by  realizing  that 
B2,  Bi,  and  Bo  are  time  varying  matrices  and  then  discretized  to  form  Equation  (2.40) 
with  a  sample  period  T: 

u(t)  =  B2(t)d{t)  +  B1(t)e(t)  +  Bo{t)6(t)  (2.39) 

u (k)  -  A2(k>T)0(k-2)  +  A1(k,T)0(k-l)  +  Ao(k,T)0(k)  (2.40) 

Assuming  a  constant  sample  period  and  using  a  discretized  joint  error  vector  definition 
of  9e(k)  —  Oj(k)  -  9(k)  allows  Equation  (2.40)  to  take  on  the  form: 

A0 (k)0e(k)  =  -u (ft)  -  Ai (k)0e(k  -  1)  -  A 2(k)ee(k  -  2)  +  A0(k)Od(k) 

+  Al(k)0d(k-l)  +  A2(k)0d(k-2)  (2.41) 
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Equation  (2.41)  suggests  that,  in  order  to  influence  the  joint  angle  error  9e  completely, 
you  need  a  control  law  of  the  general  form: 

u (*)  =  Pi(*)0e(*  -  1)  +  P 2(k)0e{k  -  2)  +  Q0 (k)9d(k)  +  Qi (k)9d(k  -  1) 

+  Q 2(k)0d(k  -  2)  +  v(k)  (2.42) 

where  Pi(fc)  and  P2(&)  are  time- varying  matrices  which  develop  the  feedback  torque  based 
upon  joint  angle  errors  and  Q0(fc),  Qi(&)>  and  Q2(fc)  are  time- varying  matrices  which, 
when  multiplied  by  the  desired  trajectory  positions,  produce  the  feedforward  torques.  The 
last  quantity,  rj(k),  is  an  auxiliary  signal  added  to  enhance  adaptation. 

Tarokh  utilized  Popov  hyperstability  to  ensure  that  the  adaptation  error  asymptot¬ 
ically  approached  zero.  A  complete  derivation  is  in  [22,  23].  Stability  conditions  were 
satisfied  by  using  proportional  plus  integral  adaptation  laws  to  derive  the  gain  matrices 
for  the  ith  joint.  The  matrices  needed  to  compute  the  control  torques  are  developed  from: 


Pii(k)  = 

Pu(k  -  1)  +  6't(k)6ei(k  -  1  )E1Pi 

+  9et(k  -  1  )Oei(k  -  2)[Eui  -  EiPi] 

(2.43) 

P2i(k)  = 

P2t{k  -  1)  +  L{k)9ci(k  -  2 )E2Pi 

+  *  i(k  -  l)9n(k  -  3 )[E2It  -  E2Pt] 

(2.44) 

Qoi(k)  = 

Qoi(k  -  1)  +  9ei(k)9di(k)FoPi 

+  9ei(k  -  l)6dt(k  -  l)[Fo/«  ~  PoPi] 

(2.45) 

Qu(k)  = 

Qu(k  -  1)  +  Oei(k)$di{k  -  1  )FlPi 

+  Oei(k  -  1  )9di{k  -  2)[F\u  -  F\P{ ] 

(2.46) 

Q2,(k)  = 

Q2i(k  -  1)  +  9C i(k)9di(k  -  2)F2Pl 

+  L(k  -  1  )0dl(k  -  3 ){F2Il  -  F2Pi] 

(2.47) 

rii(k)  = 

Vt(k  -  1)  +  PPt9ei(k )  +  {[hi  -  fipi)9ci{k  -  1) 

(2.48) 

L(k)  = 

r2iQci(k  -  1)  +  r:itOet{k ) 

r2i  - 

'■qAiiA2l(All  +  X2l) 

i':ii  -■ 

t(  1  +  Ai,A2l) 
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Table  2.2.  Tarokh’s  Design  Parameters  for  Simulations  and  Experimental  Results 


Computer  Simulation 

Experimental  Results 

Feedback  terms: 

Eip  =  01 

Ei;  =  31 

E2  p  =  01 

E2/  =  31 

Eip  =  01 
Ei;  =  21 

E2p  ~  01 
E2/  =  21 

Constants: 

c*i  =  200 

=  100 

Not  Given 

Eigenvalues: 

A2  =  -0.71 

iiMill 

A2  =  -0.51 

Auxiliary  terms: 

(3P  =  301 

f3j  =  0.021 

(3P  =  501 

(3,  =  0.51 

Initial  Conditions: 

771(0)  =  0 
Pi(0)  =  0 

772(0)  =  0 
P2(0)  =  0 

IlgjSlBm 

■aiiiBM 

772(0)  =  3 
P2(0)  =  0 

Sample  Period: 

Ts  =  1msec 

Ts  =  7msec 

where  the  E  and  F  terms  are  constant  symmetric  positive  definite  gain  matrices,  the  (3 
terms  are  positive  scalars,  and  the  subscripts  P  and  I  denote  proportional  and  integral 
terms.  The  weighted  joint  angle  errors,  9et ,  involve  constants  derived  from  positive  con¬ 
stants  aty  while  Aj,  and  A2Z  are  eigenvalues  chosen  from  the  error  reference  model  such 
that  their  magnitudes  are  less  than  1. 

2.3>4  Simulation  and  Experimental  Results  Discussion.  Like  Seraji,  Tarokh  dis¬ 
cussed  different  conditions  in  his  computer  simulation  and  experimental  results  sections. 
His  simulation  discussion  evaluated  a  3-second  trajectory  of  (90°,  0°)  to  (0°,  90°)  for  Links  2 
and  3  of  a  PUMA-560  using  only  feedback  and  auxiliary  torques.  His  control  law  was  [22]: 

u{k)  =  P i(k)6e{k  -  1)  +  P 2(k)6e(k  -  2)  +  7 j(k)  (2.49) 

This  control  law  uses  only  feedback  and  auxiliaiy  compensation.  Furthermore,  only  integral 
adaptation  was  used  to  compute  the  feedback  compensation.  The  parameters  for  the 
computer  simulation  are  shown  in  Table  2.2. 

lie  evaluated  several  other  tiajcctories,  initial  positions,  and  payload  conditions  and 
discovered  tracking  errors  of  less  than  1.0  degree.  The  simulations  also  address  the  issue 
of  adaptation  rate.  Instability  could  result  from  a  long  sample  period. 
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The  experimental  results  section  concentrated  on  three  different  test  scenarios  sum¬ 
marized  by:. 

1.  Movement  of  Joint  2  only  from  0°  to  50°  in  2  seconds  with  all  other  joints  stationary 
and  no  payload. 

2.  Movement  of  Joint  2  only  from  0°  to  50°  in  2  seconds  with  all  other  joints  stationary 
and  a  3  kg  payload. 

3.  Movement  of  Joints  2  and  3  from  (0°,0°)  to  (50°, -50°)  in  2  seconds  with  all  other 
joints  stationary  and  no  payload. 

The  parameters  for  the  experimental  runs  are  different  from  the  computer  simulations. 
They  are  also  shown  in  Table  2.2.  The  tracking  responses  of  Joint  2  comparing  test  scenario 
(1)  versus  (2)  and  scenario  (1)  versus  (3)  are  given  in  [22,  23].  These  plots  show  that  the 
different  test  conditions  cause  little  difference  in  trajectory  tracking. 

Questions  left  unanswered  by  Tarokh  which  will  be  addressed  by  this  effort  include: 

1.  Is  this  controller  applicable  under  a  variety  of  harsher  trajectory  conditions? 

2.  Are  the  parameters  chosen  for  one  trajectory  good  for  only  that  particular  trajectory 
or  several  trajectories? 

3.  Is  this  control  algorithm  sensitive  to  the  speed  of  the  trajectory? 

2Jf  Summary 

This  chapter  opened  with  a  discussion  of  the  development  of  an  AMI3C  controller  and 
detailed  the  implementation  that  was  used.  An  overview  of  decentralized  adaptive  control 
followed.  Tarokh’s  and  Seraji’s  versions  were  examined. 

The  next  chapter  will  discuss  the  results  of  adaptive  model-based  control  research.  It 
will  be  followed  by  a  discussion  of  the  decentralized  digital  adaptive  contiol  effoit^. 
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III.  Adaptive  Model-Based  Control 


3.1  Introduction 

In  several  articles  Leahy  and  Whalen  have  espoused  Adaptive  Model-Based  Control 
(AMBC)  as  a  method  for  robotic  control  [12, 13].  The  aim  of  this  thesis  effort  was  to  include 
a  third  set  of  friction  parameters  in  the  adaptive  control  algorithm.  These  parameters 
addressed  static  friction  in  an  effort  to  reduce  endpoint  errors,  especially  at  the  start  of 
the  trajectory.  Initially  the  algorithms  were  tuned  to  accept  this  change.  Several  different 
trajectories  were  traversed  under  a  variety  of  initial  conditions.  Finally  the  feasibility  of 
employing  only  the  most  critical  gravity  and  inertia  terms  along  with  the  nine  friction 
terms  in  the  AMBC  algorithm  was  explored. 

3.2  Test  Environment 

The  experimental  evaluations  performed  in  this  effort  were  conducted  under  the  AFIT 
Robotic  Control  Algorithm  Development  and  Evaluation  (ARCADE)  environment.  AR¬ 
CADE  is  hosted  on  a  VAXstation  III  and  has  both  serial  and  parallel  connections  to  the 
original  PUMA  computer  bus.  The  PUMA’s  LSI-11/73  computer  now  serves  strictly  as  a 
preprocessor.  The  angular  position  and  motor  current  information  is  passed  between  the 
preprocessor  and  the  VAXstation  via  a  DRV11-J  parallel  interface.  The  parallel  device 
driver  software  is  supplied  by  the  VAXlab  software  package.  Velocity  error  is  calculated 
by  differencing  the  desired  and  calculated  velocity  signal: 

m  =  q(k)-q{k-i)]/Ta  (3.i) 

Velocity  noise  translates  into  torque  spikes  and  is  the  major  limitation  of  algorithm  feed¬ 
back  gain,  and  therefore  tracking  performance.  A  simple  low-pass  filter  has  been  added  to 
help  reduce  these  noise  effects.  Communication  restrictions,  minimal  processing  time,  and 
nominal  clock  rate  resulted  in  a  servo  rate  of  4.5  ms  (2*22  Hz)  for  these  evaluations  [11]. 
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3.3  Algorithm  Implementation 


The  specific  structure  of  the  AMBC  algorithm  implemented  during  this  research  was 
identical  to  that  of  previous  research  on  the  PUMA.  Feedforward  torques  are  produced 
based  on  estimated  and  known  parameters. 


T  =  Tff  +  rJb 


(3.2) 


r  =  (Y(0d,6d,6d)a)  +  (Kve  +  Kpe) 


(3.3) 


Simple  PD  feedback  torque  is  used  with  the  following  gains  [12]: 


Link  # 

1 

2 

3 

Kp 

640 

1330 

360 

Kv 

72 

130 

25 

To  visualize  the  layout  of  the  links,  think  of  the  Puma  as  modeled  like  a  human  torso.  Link 
1  is  the  torso,  Link  2  is  the  upper  arm,  and  Link  3  is  the  forearm.  Joints  1-3  respectively 
are  the  waist,  shoulder,  and  elbow  [6]. 


These  torques  are  digitally  implemented  for  each  joint  by  the  following  set  of  equa¬ 
tions: 


rjb(k)  = 

Kve(k  -  1)  +  Kpe(k  -  1) 

(3.4) 

Tff(k)  = 

Yi[q<l(k),  (Jd{k),  (}d(k)]a(k)  +  Y2[q<i(k),  qd(k),  qd(k)]an(k ) 

(3.5) 

II 

t— t 

1 

<U(k  -  1)  -  [q(k  -  1)  -  q(k  -  2)] /2’s 

(3.6) 

e(k-l)  = 

5? 

i 

►—» ■ 

i 

I 

(3.7) 

where  q  and  q j  represent  the  actual  and  desired  joint  angles,  k  denotes  the  sample  instant, 
and  Ts  is  the  sample  period. 

The  regressor,  Y,  is  based  on  the  structure  of  the  manipulator  system  dynamics, 
including  reflected  actuator  inertia  and  friction  terms,  and  is  dependent  on  only  the  desired 
trajectory  positions,  velocities,  and  accelerations,  not  actual  trajectory  characteristics. 


The  columns  of  the  regressor  matrix  are  arranged  based  on  a  priori  knowledge  of 
PUMA  dynamics  [13].  The  first  three  columns  of  Y  (and  elements  of  a)  are  directly 
related  to  gravitational  torques.  Inertial  parameters  are  arranged  according  to  payload 
sensitivity  analysis.  Columns  Y(ll  -  19)  are  viscous,  Coulomb,  and  static  friction  terms. 

The  regressor  matrix  was  expanded  from  a  3  x  28  matrix  listed  in  [13]  to  a  3  X  31 
matrix  to  incorporate  the  three  static  friction  terms.  These  terms  were  inserted  as  elements 
Y(l,17)  through  Y(3,19)  as  shown: 


^ld^sign(0id) 

0 

0 

0 

e—ro&dl  sign($2d) 

0 

0 

0 

e-Mhd\ sign($3d) 

The  term  ro  is  a  lubrication  coefficient  which  is  selected  to  be  40  [5]  and  6\d  is  the  desired 
velocity  of  the  first  joint. 

The  remaining  columns  of  Y  were  shifted  to  be  Y(z,20)  through  Y(z,31)  where  i  is 
the  joint  number.  Individual  static  friction  terms  are  independent  of  other  joint  trajecto¬ 
ries,  but  other  terms  of  the  regressor  are  coupled.  Differences  in  one  link’s  performance 
attributed  to  this  compensation  may  still  indirectly  affect  the  performance  and  adaptation 
of  other  links. 

The  parameter  vector  and  regressor  matrix  are  assumed  to  have  two  components. 
The  “known”  parameters,  an,  will  remain  constant  during  program  execution  while  the 
unknown  parameters,  a,  will  be  adapted.  The  elements  of  the  regressor  matrix  correspond¬ 
ing  to  the  known  parameters  are  denoted  by  Y2.  Feedforward  torque  components  derived 
from  the  known  parameters  and  the  regressor  are  precomputed  and  stored  off-line. 

The  columns  of  the  regressor  matrix  corresponding  to  the  unknown  parameters  are 
denoted  by  Y\.  The  controller  software  is  written  such  that  the  column  size  of  Y]  is 
selected  by  the  program  operator.  The  unknown  parameters  are  initialized  to  operator- 
determined  values.  This  portion  of  the  feedforward  torque  is  calculated  on-line  [12].  "flic 


total  feedforward  torque  is  expressed  by: 


Tff  =Yi(Bd,dd,0d)a.  +  Y2(0d,0d,9d)^n  (3.8) 

The  unknown  parameters  are  adapted  by  a  Desired  Compensation  Adaptive  Law  [19]: 

a  =  J  r_1  YT(0d,  ed,  od)[(ed  -e)  +  A{ed  -  e)\  (3.9) 

This  adaptation  is  digitally  implemented  by[13]: 

a(k)=  f  r_1Y]'(gd(A:),4<i(A:),5d(L'))[e(^-l)  +  Ae(fc-l)]  (3.10) 

Jo 

The  vector  A  is  selected  to  be  the  ratio  of  the  position  feedback  loop  gain  over  the  velocity 
feedback  loop  gain.  Selection  of  T~l  is  discussed  in  a  later  section. 

3.4  Experimental  Evaluation 

3,4.1  Introduction.  Initial  testing  of  the  AMBC  algorithm  examined  the  effects  of 
adding  three  more  friction  parameters.  These  tests  were  conducted  with  an  initialized  pa¬ 
rameter  vector  and  with  an  uninitialized  vector.  The  values  used  to  initialize  the  parameter 
vector  came  from  [14].  The  parameters  were  also  rearranged  to  conduct  tests  using  only 
the  three  aforementioned  gravitational  terms,  one  inertial  parameter,  and  all  nine  friction 
terms.  Tests  were  conducted  over  a  variety  of  trajectories,  listed  in  Table  3.1,  under  zero 
payload  conditions.  Test  results  are  in  Section  3.5.  Motor  saturation  constraints  limited 
testing  with  a  2  kg  payload  to  Trajectories  1  and  3  only. 

Trajectories  1,  2,  and  3  all  experience  angular  position  movement  of  (95°,  45°,  -105°). 
Trajectory  1  undergoes  this  movement  in  1.5  seconds  instead  of  the  2.0  seconds  used  in 
Trajectories  2  and  3  and  is  the  standard  evaluation  trajectory  used  in  previous  PUMA 
evaluations  [13].  The  desired  trajectory  velocity  and  acceleration  components  will  be  iden¬ 
tical  for  Trajectories  2  and  3.  These  three  test  trajectories  will  allow  testing  to  consider 
the  effects  of  different  starting  positions  and  trajectory  speeds. 


Table  3.1.  Test  Trajectories  (degrees)[13] 


Number 

Start 

Finish 

Time  (sec) 

i 

-50,-135,135 

45,-90,30 

1.5 

2 

-50,-205,90 

45,-160,-15 

2.0 

3 

0,-180,180 

95,-135,75 

2.0 

4 

0,-90,90 

-95,-135,-15 

2.0 

5 

0,-45,135 

-95,-90,20 

2.0 

6 

0,-90,90 

95,-135,195 

1.5 

Trajectories  4  and  5  apply  identically  generated  trajectory  commands  to  different  ini¬ 
tial  conditions.  When  compared  to  Trajectories  2  and  3,  the  initial  positions  differ  and  the 
movement  of  Joints  1  and  2  is  opposite.  (Trajectories  4  and  5  move  (-95°,  -45°,  -105°)). 
The  respective  desired  position  and  acceleration  terms  should  also  differ  in  sign  from  those 
of  Trajectories  2  and  3.  These  two  test  trajectories  permit  testing  to  consider  the  effects 
of  different  starting  positions  and  direction  of  motion. 

The  movement  of  Trajectory  6  is  similar  to  Trajectory  1.  Joint  3  moves  in  the  opposite 
direction  while  Joint  1  and  2  movement  is  the  same.  Desired  velocity  and  acceleration  terms 
for  Joint  3  should  have  opposite  signs.  This  trajectory  can  be  used  to  determine  the  effects 
of  different  starting  positions,  direction  of  motion,  and  trajectory  speed. 

3.4-2  Sixteen- Parameter  Adaptation  Testing .  All  trajectories  in  Table  3.1  were  run 
to  emulate  Leahy  and  Whalen’s  adaptive  algorithm  using  16-parameters  in  the  a  vector 
[13],  validate  the  modified  software,  and  account  for  any  minor  changes  in  the  manipulator 
since  their  study.  This  was  accomplished  by  setting  the  first  three  terms  of  an,  which 
correspond  to  the  static  friction  terms,  to  zero.  Trajectories  1  and  3  were  also  run  carry¬ 
ing  a  2  kg  payload.  The  diagonal  T_1  matrix  was  set  to  Leahy  and  Whalen’s  values  of 
(120,120,120,0,90,90,90,15,150,5,80,30,15,80,80,80).  The  zero  in  the  fourth  element  denotes 
that  the  fourth  parameter  of  a  will  not  be  adapted,  but  remain  at  its  initial  value.  The 
effects  on  adaptation  caused  by  knowledge  (or  lack  thereof)  of  the  values  of  a  were  explored 
by  using  initial  approximations  of  a  and  an  equal  to  a  set  of  previously  determined  values 
(initialized)  or  a  —  an  =  0  (uninitialized). 


3.4-3  Nineteen-Parameter  Adaptation  Testing. 


3.4- 3. 1  Nineteen- Parameter  Tuning.  In  this  series  of  evaluations,  the  AMBC 
algorithm  was  used  in  conjunction  with  a  19  parameter  a  vector.  The  static  friction  terms 
were  now  part  of  the  unknown  parameter  vector.  An  a  priori  estimate  of  a(17  - 19)  =  1  was 
used  in  determining  the  best  set  of  r_1(17-  19)  for  the  adaptation  expressed  in  Equation 
(3.9).  The  first  16  values  of  T“1(l  -  16)  were  left  the  same  as  noted  in  the  last  section. 

Trajectory  1  was  run  several  times  using  different  values  for  r_1(17).  The  values 
of  r_1(18  —  19)  were  set  to  zero.  The  value  that  produced  the  best  performance  was 
selected,  then  several  more  runs  were  executed  to  determine  r_1(18),  and  finally  r_1(19). 
The  tuning  procedure  was  repeated  for  Trajectory  1  with  a  2  kg  payload.  The  values  of 
r_1(17  -  19)  were  chosen  to  be  (200,15,15). 

3. 4- 3. 2  Nineteen- Parameter  Adaptation  Runs.  All  of  the  trajectories  listed  in 
Table  3.1  were  run  using  the  T"1  matrix  obtained  from  tuning  the  AMBC  algorithm  for  19 
parameters.  The  test  conditions  included  running  the  AMBC  algorithm  with  and  without 
payloads  and: 

1.  with  a  parameters  initialized  (set  to  predetermined  values).  Then: 

(a)  Evaluate  the  algorithm  for  single-run  performance. 

(b)  Conduct  multiple  runs  with  a  set  to  values  derived  during  the  previous  runs  to 
demonstrate  the  value  of  learning  in  a  repetitive  task  environment. 

2.  with  a  parameters  uninitialized  (set  to  zero).  Then: 

(a)  Evaluate  the  algorithm  for  single-run  performance. 

(b)  Conduct  multiple  learning  runs  to  learn  the  a  parameters. 

3.44  Thirteen- Parameter  Adaptation  Testing.  The  last  set  of  tests  undertaken  us¬ 
ing  the  AMBC  algorithm  necessitated  rewriting  the  test  softwaie.  These  tests  evaluated  the 
effects  of  using  just  three  gravity  parameters,’ one  inertial  parameter,  and  all  nine  friction 


parameters  as  the  first  thirteen  elements  of  a.  The  regressor  matrix  was  again  reconfig¬ 
ured  to  place  the  friction  compensation  terms  originally  in  Y(i,  11  -  19)  into  Y(i,5  -  13), 
respectively,  where  i  is  the  joint  number.  Regressor  components  Y(i,5  -  10)  were  replaced 
as  Y(z,  14  -  19).  The  remaining  regressor  terms  were  left  in  their  initial  positions.  The 
regressor  and  manipulator  parameter  terms  were  moved  so  that  all  of  the  terms  now  as¬ 
sociated  with  the  known  parameters  could  be  set  to  their  predetermined  values  and  the 
corresponding  feedforward  torques  precalculated  and  stored. 

3. 4- 4-1  Thirteen- Parameter  Tuning.  The  AMBC  algorithm  was  now  tested 
with  a  13-parameter  a  vector.  The  T”1  parameters  pertaining  to  the  gravity  and  iner¬ 
tial  terms  were  left  at  (120,120,120,0).  These  values  were  used  by  Leahy  and  Whalen 
in  previous  tests  with  7,  10,  and  13  element  unknown  parameter  vectors  [12,  13].  As 
in  19-parameter  tuning,  an  estimate  of  a(ll  -  13)  =  1  was  used  for  the  static  friction 
terms.  Trajectory  1  was  run  several  times  using  different  values  for  r_1(5  -  7).  The 
values  of  F-1(8  -  13)  were  set  to  zero.  The  values  of  r-1(5  -  7)  that  produced  the 
best  performance  were  selected,  then  several  more  runs  were  performed  to  determine 
T“1(8  -  10),  and  finally  r““1(ll  -  13).  The  values  of  T~l(b  -  13)  were  chosen  to  be 
(600,600,200,600,600,200,200,200,200). 

3. 4- 4-2  Thirteen- Parameter  Adaptation  Runs.  All  of  the  trajectories  listed  in 
Table  3.1  and  the  two  payload  bearing  trajectories  mentioned  in  Section  3.4.1  were  run 
using  the  T”1  matrix  obtained  from  tuning  the  AMBC  algorithm  for  13  parameters.  The 
test  conditions  were  the  same  as  in  the  19-paramcter  case.  The  next  section  discusses  the 
results  obtained  from  the  16-,  19-,  and  13-parameter  tests. 

3.5  Adaptive  Model-Based  Control  Algorithm  Test  Results 

3.5.1  Introduction.  All  test  trajectories  given  in  Table  3.1  were  run  under  no-load 
conditions.  Trajectories  1  and  3  were  also  run  carrying  a  2  leg  payload.  All  test  trajectories 
were  mu  scveial  times  using  the  adaptive  parameters  derived  fiom  the  previous  learning 
run  as  initial  parameter  estimates. 
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This  discussion  compares  tests  completed  under  similar  test  scenarios.  The  13-,  16- , 
and  19-parameter  adaptive  algorithms  were  run  under  the  conditions  of: 

•  uninitialized  parameters  and  no  load 

•  uninitialized  parameters  and  a  2  kg  load 

•  initialized  parameters  and  no  load 

•  initialized  parameters  and  a  2  kg  load 

The  adaptive  algorithm  test  results  were  also  compared  to  evaluations  that  used  a  Single 
Model-Based  Controller. 

The  Single  Model-Based  Controller  (SMBC)  for  these  comparisons  was  formed  from 
the  model-based  controller  used  in  the  AMBC  algorithm.  The  SMBC  controller  used  the 
same  PD  feedback  torques  as  the  AMBC  controllers  and  feedforward  torques  based  solely 
on  the  desired  trajectory  attributes  and  assumed  manipulator  parameters  as  shown  in 
Equations  (2.3)  and  (2.20).  The  entire  manipulator  parameter  vector  a  was  initialized  to 
predetermined  values  and  the  T-1  vector  was  set  to  zero  to  prohibit  any  of  the  parameters 
from  being  adapted. 

3.5.2  Single  Run  Tests  With  Uninitialized  Parameters.  Figures  3. 1-3.3  are  represen¬ 
tative  of  the  plots  in  Appendix  D.  These  figures  show  that  the  uninitialized  13-parameter 
AMBC  algorithm,  of  which  nine  parameters  are  friction  parameters,  produces  better  tra¬ 
jectory  responses  than  the  uninitialized  16-  or  19-parameter  AMBC  controllers  for  single 
run  tests.  This  may  be  attributed  to  the  large  gains  in  the  13-parameter  T_1  vector.  The 
manipulator  parameters  may  adapt  more  quickly  than  those  determined  using  the  16-  or 
19-parameter  algorithms’  T"1  vectors.  The  16-  and  19-parameter  error  plots  trace  eacli 
other  closely  throughout  the  entire  set  of  single  run  tests  with  uninitialized  parameters. 
The  addition  of  the  three  static  friction  terms  to  the  16-parameter  AMBC  controller  does 
not  have  a  significant  influence  on  trajectory  tracking  for  single  run  uninitialized  tests. 

3.5.3  Learning  Runs  With  Uninitialized  Parameters.  The  13-,  16-,and  19-parameter 
AMBC  controllers  were  run  several  times  to  allow  the  algorithms  to  learn  the  manipulator 
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parameters.  Appendices  A,  B,  and  C  contain  plots  of  the  tracking  errors  developed  in  each 
of  the  trajectories  as  the  pc  rameters  were  being  learned.  Appendix  E  displays  comparisons 
of  the  final  results  of  learning  for  all  test  trajectories  and  provides  a  comparison  to  SMBC 
testing.  Appendices  B  and  C  plots  show  that  the  overall  tracking  errors  are  reduced  and 
peak  errors  are  reduced  by  as  much  as  80%  (see  Figures  B.10  and  C.10)  of  the  original 
errors.  Figures  3. 4-3. 6  highlight  the  results  of  learning  with  uninitialized  parameters  for 
Trajectory  4.  The  uninitialized  19-  and  16-parameter  AMBC  controllers  produce  results 
superior  to  the  SMBC  and  13-parameter  algorithms.  In  most  test  cases  the  19-parameter 
AMBC  controller  is  closely  tracked  by  the  16-parameter  controller.  In  the  remainder  of 
the  tests  the  19-parameter  controller  is  better.  The  extra  friction  terms  slightly  improve 
the  trajectory  response. 

Even  after  learning,  the  uninitialized  13-parameter  AMBC  controller  is  inferior  to  the 
SMBC  controller.  Over  half  of  the  plots  in  Appendix  A  show  that  the  learning  runs  actually 
caused  this  controller  to  produce  worse  tracking  errors  than  those  experienced  in  single  lun 
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Figure  3.4.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  2 
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Joint  3  Tracking  Errors 


Figure  3.6.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  2 
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tests.  The  13-parameter  AMBC  may  have  made  an  initial  step  in  the  wrong  direction  when 
trying  to  estimate  its  viscous  friction  parameters  and  forced  to  stay  in  that  area  because 
of  the  heavy  emphasis  on  friction  terms  imposed  by  the  corresponding  large  F-1  values. 
The  results  of  13-parameter  adaptation  with  initialized  parameters  are  much  improved  and 
the  manipulator  parameter  adaptations  are  more  in  line  with  the  initialized  values.  The 
viscous  friction  terms  are  either  positive  or  much  closer  to  zero  than  the  estimates  made  in 
the  uninitialized  cases.  Retuning  the  first  four  T*1  parameters  before  tuning  the  friction 
components  may  have  relieved  some  of  this  error,  and  tuning  during  learning  runs  instead 
of  single  trajectory  testing  could  help  improve  this  algorithm’s  performance. 

Appendices  L  through  N  contain  tables  listing  the  parameters  learned  through  repet¬ 
itive  testing..  Appendix  0  provides  a  side-by-side  comparison  of  comparable  parameters. 
Table  3.2  shows  the  results  of  no-load  Trajectory  1  testing.  The  three  columns  on  the 
right  contain  the  manipulator  parameters  deteimined  after  the  ninth  learning  runs  of  the 
uninitialized  parameter  tests.  These  algorithms  should  be  trying  to  find  the  same  a  vahrs. 
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Table  3.2.  Learned  Parameter  Values  -  Trajectory  1 


Trajectory  1  Initialized 

Trajectory  1  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

i 

0.30 

0.85 

1.04 

0.87 

-4.96 

5.81 

6.78 

2 

-52.09 

-54.86 

-53.83 

-53.15 

-30.73 

-17.64 

-17.88 

3 

-7.53 

-7.30 

-8.81 

-9.63 

-10.12 

-8.46 

-8.88 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.23 

0.38 

-1.08 

-1.18 

6 

-1.57 

-2.86 

-2.77 

-3.34 

-3.36 

7 

2.10 

2.95 

3.40 

0.65 

0.82 

8 

-0.03 

-0.70 

-0.65 

-1.20 

-1.24 

9 

0.67 

2.99 

3.04 

11.41 

11.58 

10 

-0.12 

-0.23 

-0.16 

-1.90 

-1.80 

11 

4.50 

-1.73 

4.81 

4.23 

-36.53 

7.04 

7.00 

12 

3.50 

-0.23 

4.08 

4.16 

-23.76 

3.58 

13 

3.50 

1.09 

3.45 

3.52 

-9.25 

3.18 

14 

5.95 

7.62 

7.46 

7.35 

10.42 

5.02 

4.74 

15 

6.82 

10.58 

9.51 

9.88 

20.56 

17.00 

17.13 

16 

3.91 

4.30 

5.86 

6.02 

4.72 

9.45 

17 

1.00 

0.21 

1.64 

-5.11 

1.22 

18 

1.00 

-0.34 

0.87 

-2.43 

-0.16 

19 

1.00 

0.37 

0.91 

-5-90 

-0.33 

The  manipulator  parameters  found  by  the  16-  and  19-parameter  AMBC  controllers  are 
similar  from  top  to  bottom.  There  is  no  correlation  between  these  manipulator  parame¬ 
ters  and  those  developed  by  the  13-parameter  AMBC  controller.  The  discrepancies  of  the 
13-parameter  controller  may  have  been  caused  by  using  fewer  parameters  to  produce  the 
same  level  of  robot  control  and  needing  different  values  to  attain  the  goal,  or  by  the  large 
T-1  vector  causing  the  manipulator  parameters  values  to  change  much  more  drastically 
than  in  the  16-  or  19-parameter  cases.  The  13-parameter  algorithm  suffers  severely  when 
half  of  its  parameters  are  given  an  incorrect  sign  as  a  result  of  adaptation.  Other  tables  in 
Appendix  0  show  than  the  uninitialized  19-  and  16-parameter  controllers  produce  results 
similar  to  each  other  for  each  trajectory.  The  parameter  values  change  from  trajectory  to 
trajectory  which  is  consistent  with  the  findings  of  Leahy  and  Whalen  [13]. 
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3.5.4  Single  Run  Tests  With  Initialized  Parameters .  Trajectory  errors  resulting 
from  testing  with  initialized  parameters  are  at  least  50%  smaller  than  those  experienced 
running  single  run  tests  with  uninitialized  parameters.  Figures  3. 7-3. 9  are  representative 
of  the  no-load  trajectories  in  Appendix  I.  For  every  no-load  trajectory  test  run,  the 
13-parameter  AMBC  algorithm  produced  the  best  response,  except  Trajectory  6  -  Joint  3 
where  large  oscillations  developed  in  the  middle  of  the  trajectory  tracking  error.  The  16- 
and  19-parameter  models  develop  tracking  responses  that  are  similar  to  each  other.  As 
expected,  the  SMBC  controller  was  ineffective  when  compared  to  any  of  these  AMBC 
controllers. 

The  2  kg  payload  test  runs  all  experienced  oscillations  and  degradation  over  the  no- 
load  tests.  The  16-parameter  algorithm  developed  the  most  favorable  set  of  error  profiles. 
The  19-parameter  system  was  a  very  close  second  and  the  13-parameter  model  third.  The 
13-parameter  algorithm  may  have  been  penalized  for  its  smaller  set  of  adapting  parame¬ 
ters.  The  single  run  tests  utilizing  initialized  parameters  were  nearly  comparable  to  the 
uninitialized  parameter  cases  after  learning.  These  sets  of  initialized  runs  show  that  there 
is  no  substitute  for  good  initial  estimates. 

3.5.5  Learning  Runs  With  Initialized  Parameters.  AMBC  controllers  with  13,  16, 
and  19  initialized  parameters  were  allowed  to  use  the  parameters  learned  during  one  run  as 
the  initial  parameter  estimate  for  the  next  run.  Appendices  F,  G,  and  H  contain  plots  of  the 
errors  generated  during  the  learning  process.  As  was  the  case  with  uninitialized  parameters, 
learning  runs  for  the  initialized  13-parameter  AMBC  controller  did  not  necessarily  mean 
that  the  trajectory  tracking  response  was  going  to  improve.  Several  of  the  figures  in 
Appendix  F  show  that  the  errors  increased  as  the  controller  was  supposed  to  be  learning 
the  manipulator  parameters.  Neither  the  16-  nor  19-paramcter  AMBC  controller  could 
improve  tracking  response  for  Joint  3  of  Tiajectory  6.  These  controllers  did  improve  or 
maintain  the  tracking  responses  for  all  of  the  other  test  cases.  Some  of  the  tracking  error 
plots  in  Appendices  G  and  II  showed  peak  error  reductions  of  more  than  60%. 

The  learning  process  was  complete  by  the  fifth  run  for  ail  three  algorithms.  Figures 
3.10-3.12  show  one  set  of  plots  contained  in  Appendix  J.  The  appearance  of  oscillations  in 
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some  plots  is  due  to  the  small  range  covered  by  the  vertical  axis  of  the  plots.  Figure  3.11 
covers  an  entire  range  of  .003  radians  from  top  to  bottom. 

The  16-  and  19-parameter  algorithms  usually  trace  similar  error  profiles.  The  initial¬ 
ized  13-parameter  AMBC  controller  is  surprisingly  capable.  It  provided  the  best  control 
for  all  of  the  joints  of  the  robot  for  Trajectory  3..  This  controller,  however,  was  the  worst 
controller  for  the  payload  bearing  and  high  speed  trajectories.  The  19-parameter  controller 
worked  best,  and  outperformed  the  16-parameter  controller,  cn  the  heavier  Links  1  and  2. 
The  16-parameter  controller  was  best  utilized  for  Joint  3. 

The  three  columns  in  the  center  of  Table  3.2  contain  the  results  of  no-load  Trajectory 
1  testing  for  learning  with  initialized  parameters.  Again  the  manipulator  parameters  found 
by  the  16-  and  19-parameter  AMBC  controllers  are  similar.  There  are  sign  changes  in  some 
of  tlm  parameters  developed  by  the  13-parameter  AMBC  controller.  It  has  decided  that  the 
viscous  friction  terms  for  Joints  1  and  2  should  be  negative.  Tables  0.4,  0.5,  and  0.8  show 
the  16-  and  19-parameter  controllers  determining  that  parameter#!  as  a  negative  number. 
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Figure  3.10.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  4 


Joint  2  Tracking  Errors 


Time  (seconds) 


Figure  3.1 


Comparison  of  Initializer 


IB53 

19  Parameters 

13  Parameters 

■991 

16  Parameters 

Runs  After  Learning  -  Trajectory  4 


3-17 


Joint  3  Tracking  Errors 


Figure  3.12.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  4 
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The  nominal  values  were  positive.  Similar  changes  in  the  signs  of  several  parameters  are 
probably  responsible  for  the  reduction  in  the  response  for  the  13-parameter  AMBC. 

3.5.6  Comparison  of  Initialized  and  Uninitialized  Learning  Runs .  Figures  3.13-3.15 
from  Appendix  K  illustrate  the  value  of  initialized  parameters  over  uninitialized  parameters 
for  16-  and  19-parameter  controllers.  The  13-parameter  controller  is  inferior  to  these  two 
because  of  its  shortcomings  for  the  high-speed  and  loaded  tests.  These  figures  show  that  the 
tracking  errors  can  be  reduced  to  approximately  .001  radian.  If  computational  loading  is  an 
issue,  the  16-parameter  controller  performs  nearly  as  well  as  the  19-parameter  controller. 
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Figure  3.13.  Compariscn  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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Figure  3.14.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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Figure  3 .15.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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3.6  Summary 

This  chapter  discussed  the  evaluation  of  an  Adaptive  Model-Based  Controller.  Tests 
were  conducted  over  a  variety  of  trajectories  with  different  numbers  of  initialized  and 
uninitialized  parameters.  Uninitialized  tests  evaluated  the  algorithms’  ability  to  develop 
manipulator  parameters  in  the  absence  of  any  prior  information  about  the  parameters. 
The  addition  of  extra  friction  parameters  did  not  improve  the  response  of  the  basic  16- 
parameter  AMBC  during  single  run  testing.  Learning  tests  did  reduce  the  peak  errors  of 
the  single-run  tests  by  as  much  as  80%. 

Initialized  tests  provided  the  controller  with  an  a  priori  estimate  of  manipulator  pa¬ 
rameters.  The  controllers  were  also  able  to  compute  an  additional  feedforward  torque  com¬ 
ponent  based  on  known  non-adapting  parameters.  The  parameter  knowledge  helped  reduce 
tracking  errors  by  50%  with  respect  to  comparable  uninitialized  tests.  The  13-parametei 
AMBC  emphasizes  friction  components.  It  produced  the  best  single-run  responses,  but 
was  penalized  for  its  smaller  parameter  size  and  nearly  diagonal  structure  of  the  regres- 
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sor  matrix  in  the  learning  tests.  The  16-  and  19-parameter  AMBCs  reduced  most  peak 
tracking  errors  to  less  than  one  milliradian.  There  was  no  decided  advantage  to  using  19 
parameters  over  16  for  initialized  testing,  and  therefore  no  need  to  implement  the  static 
friction  parameters, 
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IV.  Evaluation  of  Decentralized  Digital  Control  Algorithms 


4 A  Introduction  to  the  Discussion  of  Tarokh’s  Algorithm 

In  a  pair  of  publications,  Tarokh  has  presented  the  development  of  a  Popov 
hyperstability-based  discrete-time  adaptive  control  scheme  [22,  23].  These  papers  were 
discussed  in  Chapter  2.  His  initial  evaluations  were  conducted  over  trajectories  that  moved 
Joints  2  and  3  of  a  PUMA-560  manipulator  from  (0°,0°)  to  (50°,  -50°).  The  time  of  the 
trajectory  was  2  seconds.  Attempts  were  made  to  validate  those  results.  Restrictions  im¬ 
posed  on  our  Puma  led  to  our  test  trajectories  starting  at  (-180°,  180°)  and  finishing  at 
(-130°,  130°).  The  design  parameters  Tarokh  selected  are  listed  in  Table  4.1. 


Table  4.1.  Tarokh’s  Design  Parameters 


Feedback  terms: 

Eip  =  0 

N5 

II 

O 

El/  =  212x2 

E2/  =  2I2x2 

Fqp  =  0 

Feedforward  terms: 

Fip  =  0 

HBfl 

F2P  =  0 

1B1 

Constants: 

(*2  =  200 

a3  =  100 

Eigenvalues: 

Ai  =  -0.512x2 

A2  =  —0.512x2 

Auxiliary  terms: 

Pp  =  5012x2 

/3;  =  0.512x2 

Initial  Conditions: 

772(0)  =  -25 

773(0)  =  3 

The  subscripts  for  the  r/  and  a  terms  refer  to  the  joint  number.  The  subscripts  on  the 
E  and  F  matrices  denote  which  P  or  Q  matrix  of  Equation  (2.42)  they  are  used  to  create. 
The  implementation  of  Equation  (2.42)  is  discussed  in  the  next  section.  The  values  of  the 
constant  a  terms  were  not  given  in  the  section  of  Tarokh’s  paper  covering  experimental 
results  [22].  The  a  values  listed  in  Table  4.1  are  taken  from  his  computer  simulation 
section.  All  values  in  this  table  are  to  be  specified  by  the  designer  to  meet  the  needs  of 
the  particular  manipulator.  In  the  attempt  to  reproduce  Tarokh’s  plots,  the  a  values  used 
were  (100,50)  and  the  A  values  were  changed  to  -0.8.  The  initial  conditions  listed  in  Table 
4.1  were  used.  If  this  control  algorithm  is  not  dependent  on  its  trajectory,  initial  conditions 
should  apply  to  any  trajectory  selected.  Feedforward  compensation  terms  were  not  listed 
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in  either  the  computer  simulation  or  experimental  results  section  and  were  therefore  taken 
to  be  zero. 

The  evaluation  runs  were  accomplished  under  the  same  conditions  for  adjacent  link 
movement  and  payload  parameters  as  Tarokh  used.  Figure  4.1  shows  the  desired  and 
actual  trajectories  covered  by  the  Joint  2  of  the  manipulator  under  the  conditions  of  no 
movement  in  Joint  3,  with  and  without  payload.  Figure  4.2  shows  the  desired  and  actual 
Joint  2  trajectories  with  and  without  Joint  3  movement  under  no-load  condition.  The  plots 
of  actual  trajectory  errors  shown  in  Figures  4.3  and  4.4  are  more  revealing.  These  plots 
show  how  large  the  tracking  errors  produced  really  are. 

Although  the  trajectories  generated  by  our  manipulator  do  not  demonstrate  the  same 
smoothness  shown  by  Tarokh’s  experiments,  Figures  8  and  9  of  [22],  they  supplied  enough 
incentive  to  explore  the  algorithm  further  in  the  hope  of  attaining  suitable  parameters  to 
provide  accurate  trajectory  tracking.  The  next  section  describes  the  approach  taken  in 
tuning  and  testing  Tarokh’s  algorithm. 

4-2  Examination  of  Tarokh’s  Algorithm 

The  discrete-time  adaptive  control  algorithm  implemented  in  this  effort  is  given  in 
Equation  (2.42)  as: 

u(k)  =  Pi(k)Qe(k  -  1)  +  J>2(k)0e(k  —  2)  +  Qo(k)0d(k)  +  Q\(k)Od(k  —  1) 

+  W)0d{k-2)  +  ri(k)  (4.1) 

The  control  algorithm  can  be  divided  into  feedback,  feedforward,  and  auxiliary  torque 
components  as  follows: 

Feedback  Torque:  Pi(fc)0e(Jfe  -  1)  +  P2(k)Oe(k  -  2) 

Feedforward  Torque:  Qo (k)Qd(k)  +  Qi (k)Od(k  —  1)  +  Q2 (k)&d(k  -  2) 

Auxiliary  Torque:  77 ( / c ) 
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Joint  2  Trajectories 


Figure  4.1.  Joint  2  Trajectories  with  no  Joint  3  Movement 
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figure  4.2.  Joint  2  Trajectories  with  no  Payload 
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Figure  4.3.  Joint  2  Errors  with  no  Joint  3  Movement 
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where  the  components  are  modeled  by: 


P«(A) 

=  P x(k  -  1)  +  9e(k')9e(k  -  x)~Exp 

+  9e(k  -  1  )9e(k  -  x  -  1)[E,/  -  E*p]  x  =  1,2 

(4.2) 

QyW 

=  Qy(^  -  1)  +  9e(k)9d(k  -  y)Fyp 

+  9e(k  -  1  )9d(k  -y-  l)[Fy/  -  FyP]  y  =  0,1,2 

(4.3) 

v(k) 

=  ri(k  - 1)  +  I3p9e(k)  +  (f3j  -  /3p)9f(k  -  1) 

(4.4) 

The  weighted  joint  angle  error  6e  is  formed  from  the  matrix  parameters  a,  A1?  and  A2  by: 


R-2 

—  a:AiA2(Ai  +  A2) 

(4.5) 

r3 

=  <*(^3x3  +  A*  A2) 

(4.6) 

9e(k) 

=  R 20e(*  -  1)  +  R3 0e(k) 

(4.7) 

The  terms  a,  A,  R2,  and  R3  are  constant  diagonal  matrices  and  the  weighted  errors  are 
common  to  all  torque  components. 

The  torque  components  were  ranked  in  descending  order  of  magnitude,  assuming  the 
presence  of  small  joint  angle  errors,  to  determine  where  to  start  the  process  of  developing 
parameter  values  for  our  particular  controller.  The  ordei  and  reasoning  are: 

1.  Auxiliary  torque  -  This  term  is  the  weighted  sum  of  the  two  most  recent  weighted 
errors  and  is  only  first  order  as  far  as  the  error  terms  are  concerned. 

2.  Feedforward  torques  -  The  three  feedforward  torque  components  are  products  of  first 
order  error  terms  and  second  order  angular  position  terms.  Squaring  these  quantities 
eliminates  dependence  on  the  sign  of  the  position  terms.  The  position  errors  are 
assumed  to  be  small  while  the  position  terms  can  be  as  large  as  ±tt.  This  term  could 
cause  problems  if  its  components  grow  too  large. 

3.  Feedback  torques  -  These  two  torques  arc  third  order  products  of  the  position  errors. 
If  errors  are  in  the  range  of  0.005  radians,  these  products  are  on  the  order  of  10“'. 


For  the  desired  small  errors,  the  parameters  associated  with  these  terms  may  need  to 

be  extremely  large  to  produce  any  noticeable  torques. 

4*3  The  Tuning  Process 

Tarokh  stated  that  this  algorithm  includes  several  parameters  to  be  chosen  by  the 
designer.  In  this  case  there  are  design  parameters  and  6  initial  conditions  to  be  selected 
for  each  of  three  links.  In  view  of  the  restrictions  imposed  by  a  decentralized  controller, 
all  design  parameters  are  viewed  as  vectors  or  diagonal  matrices.  This  algorithm  is  not 
based  on  the  manipulator’s  dynamics  model,  as  are  adaptive  model-based  algorithms.  The 
parameters  are  selected  solely  on  the  basis  of  trial  and  error. 

4.3.1  Weighted  Error  Terms.  The  first  parameters  selected  for  each  link  were  Ai  and 
A2  over  a  2-second  trajectory  of, (-50°,  -135°,  135°)  to  (45°,  -90°,  30°).  This  trajectory  is 
denoted  as  Trajectory  0.  The  first  runs  were  conducted  with  the  feedforward  component 
of  a  Single  Model-Based  Controller  (SMBC)  added  to  limit  the  errors  experienced  during 
the  heuristic  process  of  determining  an  initial  set  of  A  parameters.  The  a  and  j3  terms 
used  in  veiifying  the  algorithm  >vere  chosen  until  better  values  could  be  identified.  Values 
of  Ai  and  A2  were  tested  for  each  joint  in  the  manner  of: 

(ai,A2)  =  ±(C.2i+0.1,0.2j  +  0.i)  f  =  0,1, 2, 3, 4  j  =  0,1, 2, 3,4 

to  test  values  all  over  che  range  o'f  |A|  <  1.  Points  of  (Ax,  A2)  in  the  fourth  quadrant  form  the 
same  weighted  errors  as  the  corresponding  points  in  the  second  quadrant.  Therefore,  the 
fourth  quadrant  was  r.ot  investigated.  The  third  quadrant  (both  Ax  and  A2  <  0)  provided 
the  best  results.  Small  magnitudes  of  A  produced  oscillations  in  the  manipulator  trajectory 
while  larger  magnitudes  iesulted  in  larger  eriors  over  smoother  trajectories.  Values  of  A 
that  produced  a  compromise  of  these  opposing  conditions  were  selected.  This  hinted  at 
the  need  for  readjusting  A  values  as  parameters  were  added. 

The  constant  a  terms  were  recognized  as  being  only  multipliers  that  increase  the 
magnitude  of  the  R2  and  R;$  terms  used  in  calculating  the  weighted  errors.  A  check  of 
Equations  (4. 2-4. 4)  reveals  that  each  P,  Q,  and  rj  term  is  a  first  order  function  with  respect 


to  de(k)  or  de(k  -  1).  All  a  values  were  set  to  1.  The  constants  for  /3,  E,  and  F  terms 
will  all  be  larger  (a  times  larger)  to  compensate  for  this. 

4.3.2  Auxiliary  Torque  Terms.  The  auxiliary  torque  component  was  identified  as 
the  largest  contributor  to  total  torque,  therefore,  the  (3P  and  f3j  terms  were  selected  first. 
These  terms  seemed  to  influence  opposite  ends  of  the  trajectory  for  their  particular  link. 
In  the  case  of  Links  1  and  2,  increasing  (3P  had  the  effect  of  reducing  the  initial  peak  errors 
while  inducing  oscillations  in  the  final  part  of  the  trajectory.  Increasing  the  integral  term, 
/3;,  quelled  these  oscillations  but  reinforced  the  initial  peak  errors.  Link  3  demonstrated 
the  opposite  behaviors  for  these  terms.  These  tradeoffs  were  considered  in  selecting  f3P 
and  (3j. 

Initial  conditions  for  the  auxiliary  input  were  introduced  to  examine  their  effect.  The 
SMBC  feedforward  term  was  removed  during  the  (3  tuning  procedure  as  soon  as  a  set  of 
parameters  produced  a  reasonable  response.  Figures  4.5  -  4.7  show  the  trajectory  errors 
incurred  during  the  steps  of  this  tuning  process.  The  oscillatory  nature  of  the  tracking 
error  plots  is  due  to  the  lack  of  velocity  gains  in  Tarokh’s  algorithm.  Without  velocity 
gain,  the  algorithm  is  suitable  for  regulation,  not  tracking  [10].  After  the  best  set  of  (3 
terms  was  selected,  the  next  step  was  to  go  back  and  adjust  the  A  terms  for  smoother 
tracking  and  then  recheck  the  (3  terms.  The  parameters  selected  during  each  step  of  the 
tuning  process  are  listed  in  Table  4.2. 

4.3.3  Feedforward  Torque  Terms.  Feedforward  torque  is  the  next  largest  component 
in  the  total  control  torque.  It  consists  of  three  components  based  on  weighted  errors 
and  angular  positions.  The  error  and  position  terms  have  time  delays  ranging  from  0 
to  3  sample  periods.  All  three  proportional  feedforward  torque  parameters,  Fop,  Fip, 
and  F2P,  were  selected  before  the  integral  parameters.  Again,  the  tradeofF  in  selecting 
sets  of  these  parameters  was  the  balance  between  reduced  initial  errors  versus  endpoint 
oscillations.  The  first  two  parameters  produced  improvements  in  manipulator  tracking. 
The  third  term,  F2P,  only  acted  to  increase  oscillations.  It  was  set  to  zero.  The  A  and  (3 
terms  were  re-adjusted  and  followed  by  a  re-evaluation  of  Fqp  and  F\p. 
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Figure  4.5.  Error  Plots  During  Tuning  Tarokh’s  Algorithm 
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Figure  4.6.  Error  Plots  During  Tuning  Tarokh’s  Algorithm 
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Figure  4.7.  Error  Plots  During  Tuning  Tarokh’s  Algorithm 
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Next,  the  integral  feedforward  parameters,  Fo /,  Fi/,  and  F2/,  were  addressed.  As  in 
the  case  of  the  proportional  parameters,  only  the  first  two  terms  were  effective.  The  third 
was  set  to  zero.  Initial  conditions  for  the  three  feedforward  components  did  not  improve 
performance.  All  previously  selected  parameters  were  re-adjusted  for  optimum  response. 

4.3.4  Feedback  Torque  Terms.  Although  these  terms  (E xp,  E uy  E 2p.  and  E 2/)  are 
included  in  the  derivation  and  demonstration  of  this  algorithm,  they  did  not  influence 
the  tracking  and  endpoint  accuracy  even  when  set  as  high  as  106.  These  terms  and  the 
feedback  torque  initial  conditions  were  set  to  zero. 

4.4  Examination  of  Seraji’s  Algorithm 

This  section  is  a  cursory  glance  at  the  Seraji's  Lyapunov-based  decentralized  adaptive 
control  algorithm  [20]  that  was  discussed  in  Chapter  2.  Leahy  also  experimented  with 
this  algorithm  [10]  and  expressed  this  controller  as  being  divided  into  feedforward  (t//), 
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Table  4.2.  Design  Parameters  Selected  During  Tuning  Tarokh’s  Algorithm 


Parameters  Under  Consideration:' 

(3p  and  (3j 

Fop  and  Fip 

Fo /  and  Fj/ 

Constants: 

a 

diag[l,l,l] 

diag[l,l,l] 

diag[l,l,l] 

Auxiliary  Terms: 

fit 

diag[10,10,2] 

diag[5,5,10] 

diag[5,5,10] 

Pp 

diag[1000, 3000, 400] 

diag[1600, 4500, 700] 

diag[1600, 4500, 750] 

Eigenvalues: 

Ax 

diag[-7.,-.7,-.7] 

diag[-.85,-.85,-.85] 

diag[-.85,-.85,-.85] 

a2 

diag[-.95,-.95,-.95] 

diag[-.9,-.9,-.9] 

diag[-.95,-.95,-.95] 

Feedforward  Terms: 

Fof 

diag[0,0,0] 

diag[1200, 1000, 500] 

diag[1050,850,525] 

Fip 

diag[0,0,0] 

diag[900,700,450] 

diag[1400, 1200, 900] 

For 

diag[0,0,0] 

diag[0,0,0] 

diag[ll,8,9] 

Fi/ 

diag[0,0,0] 

diag[0,0,0] 

diag[5,5,5] 

F2 p  and  F2 / 

diag[0,0,0] 

diag[0,0,0] 

diag[0,0,0] 

Feedback  Terms: 

Eip  and  E2p 

diag[0,0,0] 

diag[0,0,0] 

diag[0,0,0] 

En  and  E2/ 

diag[0,0,0] 

diag[0,0,0] 

diag[0,0,C] 

Initial  Conditions: 

77(0) 

diag[10,30,-7] 

diag[10,30,-7] 

diag[10,30,-7] 

feedback  (r/j),  and  auxiliary  (rax)  torque  components. 

r  =  Tnx  +  Tff  +  Tjb  (4.8) 

where: 

Tff  =  COri  +  B  -f  Add  (4.9) 

r  jb  =  Kve  +  Kpe  (4.10) 

6d  is  the  desired  trajectory  and  e  -  (9j  -  6).  A,  B,  C,  K„,  and  Kp  are  all  nxn  adaptive 
matrices  [10].  These  equations  arc  digitally  implemented  by: 


7*ar(^)  “ 

rax(k  -  1)  +  a(l)^[R(fc)  +  R(k  -  1)] 

+  f3(l)[R(k)  -  R(k  -  1)] 

(4.11) 

K r(k)  = 

Kp(k  -  1)  +  a(2)^[R (k)e(k)  +  R(k  -  1  )e(k  -  1)] 

+  /3(2)[R (k)e(k)  -  R(k  -  1  )e(k  -  1)] 

(4.12) 

K  „(*)  = 

Kv(k  -  1)  +  a(3)^[R(k)e(k)  +  R (k  -  l)e(k  -  1)] 
z 

+  /3(3)[R (k)e(k)  -  R (k  -  1  )e(k  -  1)] 

(4.13) 

A(k)  = 

A (k  -  1)  +  a(4)^[R (k)0d(k)  +  R (k  -  1  )0d(k  -  1)] 

+  j3(4)[R (k)6d(k)  -  R(k  -  1  )6d(k  -  1)] 

(4.14) 

B(t)  = 

B(k  -  1)  +  a(5)^[R(fc)0d(A:)  +  R(k  -  1  )Gd(k  -  1)] 

+  (3(b)[R{k)9d{k)  -  R (k  -  l)Od(k  -  1)] 

(4.15) 

C(k)  = 

C (k  -  1)  +  a{6)j[R(k)ed(k)  +  R(k  -  1  )6d(k  -  1)] 

+  /3(6)[R(k)6d(k)  -  R{k  -  l)Od{k  -  1)] 

(4.16) 

R(/:j  = 

Wpe(k)  +  Wve{k) 

(4.17) 

Since  this  algorithm  is  defined  as  a  decentralized  controller,  all  matrices  are  forced  to  be 
diagonal.  The  diagonal  6}  a,  and  7  matrices  used  in  Equations  (2.34-2.36)  were  imple¬ 
mented  as  constants  and  gathered  into  a  new  ol  vector.  Similarly,  the  corresponding  p, 
/3,  and  A  matrices  are  now  expressed  as  a  new  (3  vector.  Tnese  a  and  (3  vectors  and  the 
diagonal  matrices  Wp  and  Wv  are  operator-specified  constants. 

Again,  as  with  Tarokh’s  algorithm,  the  torque  components  were  ranked  in  descending 
order  of  magnitude,  assuming  the  presence  of  small  joint  angle  and  velocity  errors,  to 
determine  which  parameters  should  contribute  the  most  to  developing  the  torques  for  our 
particular  controller.  The  order  of  these  components  is  the  same  as  the  corresponding 
portions  of  Tarokh’s  controller: 

1.  Auxiliary  torque  -  This  term  is  the  weighted  sum  (and  difference)  of  the  two  most 
recent  weighted  errors  R  and  is  only  first  order  as  far  as  the  error  terms  are  concerned. 
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2.  Feedforward  torques  -  The  three  feedforward  torque  components  are  products  of  first 
order  error  terms  and  second  order  angular  position,  velocity,  and  acceleration  terms. 
The  position  errors  are  assumed  to  be  small  while  the  position  terms  can  be  as  large 
as  ±7 r.  High  instantaneous  acceleration  and  velocity  terms  must  also  be  regarded 
with  caution.  This  may  be  why  Seraji  chose  to  omit  the  feedforward  terms  from  his 
implementation  of  the  algorithm  in  [20]. 

3.  Feedback  torques  -  As  in  Tarokh’s  algorithm,  these  torques  are  third  order  products 
of  angular  property  errors.  For  the  desired  small  errors,  the  multiplicative  products 
of  the  parameters  associated  with  these  terms  may  need  to  be  extremely  large  to 
produce  any  noticeable  contribution. 

Jf./f.  1  Differences  in  Tarokh’s  and  Scraji’s  Algorithms.  There  are  two  major  differ¬ 
ences  in  these  two  decentralized  adaptive  controllers  other  than  the  stability  theorem  used. 
One  has  to  do  with  the  trajectory  components  used.  The  other  concerns  the  nature  of  the 
adaptive  matrices. 

Tarokh’s  controller,  Equation  (4.1),  is  formulated  on  only  angular  position  and  po¬ 
sition  error  terms.  All  of  the  torque  components  and  adaptive  parameters  are  based  on 
functions  of  position  and  position  error.  Scraji’s  controller,  Equations  (4.8-4.10),  uses  er¬ 
ror  terms  and  desired  trajectory  attributes.  The  adaptive  matrices  depend  on  position 
and  velocity  errors  contained  in  the  weighting  matrix  R.  These  velocity  errors,  which  are 
calculated  by  e(k )  =  (c(&)  -  e(k  -  1  )]/Ts,  where  e(k)  =  Qd(k)  -  0(k ),  are  essentially  a 
weighted  position  error  term.  The  feedforward  torque  components  and  their  respective 
adaptive  matrices  rely  on  desired  trajectory  components  which  avoids  the  use  of  noisy 
acceleration  measurements. 

In  Seraji’s  algorithm,  as  implemented  by  both  Seraji  [20]  and  Leahy  [10],  only  the 
weighted  error  terms  allow  the  matrices  to  adapt  on  a  joint  by  joint  level.  The  a  and 
(3  terms  are  constants  applied  to  a  particular  matrix.  The-  adaptive  elements  of  Tarokh’s 
controller  are  allowed  the  additional  freedom  of  diagonal  matiices  E  and  F  to  aid  in  the 
adaptation  of  the  matrices  used  in  the  feedforward  and  feedback  torque  terms.  Both 
controllers  express  the  auxiliary  term  as  a  linear  function  of  weighted  on  or  terms. 
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4.4-2  Experimental  Runs  Made  Using  Seraji’s  Algorithm.  The  set  of  tuning  param¬ 
eters  used  by  both  Seraji  [20]  and  Leahy  [10]  is  given  in  Table  4.3.  All  initial  conditions 
were  set  to  zero.  These  parameters  were  used  to  run  all  proposed  test  trajectories  listed  in 
Table  3.1.  These  parameters  reduced  the  controller  to  using  only  auxiliary  and  feedback 
torques.  Setting  a(3)  —  a(6)  =  0  deleted  the  feedforward  torque  terms.  Setting  the  /3 
terms  all  to  zero  removed  the  contribution  of  the  terms  associated  with  the  differences  of 
the  two  most  recent  weighted  error  terms.  Neither  Seraji  nor  Leahy  tuned  the  controller 
for  optimum  response.  None  of  the  experimental  results  are  particularly  impressive.  Nor 
is  it  surprising  that  the  plots  contained  in  Appendix  Q  illustrate  that  this  controller  is 
sensitive  to  initial  positions,  the  signs  of  desired  trajectory  parameters,  and  payload. 


Table  4.3.  Seraji’s  Design  Parameters 


Auxiliary  terms: 

o(l)  =  50 

/3(1)  =  0 

Feedback  terms: 

it  it 
"oo  Is? 

II  II 

00  ►— 
o  o 
o  o 

m  =  0 

m  =  o 

Feedforward  terms: 

m  =  o 

13(5)  =  0 

m  =  o 

Weighted  Error  terms: 

Wp{  1)  =  30 
Wp(  2)  =  40 
WP(3)  =  12 

Wv(l)  =  20 

Wv(2)  =  20 
Wv(  3)  =  4 

4 ‘4.3  Evaluation  of  Tarokh  and  Seraji  Controllers.  The  decentralized  adaptive  con¬ 
trollers  will  be  evaluated  on  their  ability  to  respond  to  conditions  other  than  their  tuning 
trajectory  conditions.  Is  this  controller  sensitive  to  payload,  trajectory  speed,  or  different 
trajectories? 

Tarokh’s  algorithm  was  tested  over  the  trajectories  listed  in  Table  3.1.  The  controller 
was  tuned  for  Trajectory  0.  The  test  trajectories  demonstrate  the  response  of  the  manipu¬ 
lator  to  scenarios  other  than  the  tuning  conditions,  f  igures  1.8-4.10  show  how  the  Taiokh 
and  Seraji  decentralized  adaptive  control  schemes  match  up  against  each  other.  Appendix 
Q  contains  the  set  of  plots  comparing  all  trajectories.  Appendix  P  contains  plots  which 


4-14 


Joint  3  Tracking  Errors 
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Figure  4.10.  Comparison  of  Tarokh’s  and  Seraji’s  Control  Algorithms  -  Trajectory  1 


Tarokh  Controller 

ms 

Seraji  Controller 

compare  the  error  plots  generated  from  the  tuning  Trajectory  0  to  those  for  the  other 
test  trajectories  using  the  same  design  parameters.  Trajectories  0,  2,  and  3  all  experience 
movement  of  (95°,45°,  —105°)  in  2.0  seconds  and  are  subject  to  identical  desired  velocity 
and  acceleration  commands.  Tracking  errors  for  Trajectories  2  and  3  are  larger  than  those 
for  Trajectory  0.  Joints  2  and  3  (especially  3)  go  into  oscillation  during  Trajectory  2. 
Figures  4.11-4.13  illustrate  the  Tarokh  controller’s  sensitivity  to  position  terms. 

Trajectories  4  and  5  apply  identically  generated  trajectory  commands  to  different 
initial  conditions.  When  compared  to  Trajectory  0,  the  initial  positions  differ  and  the 
desired  velocity  and  acceleration  terms  should  be  opposite  in  sign  since  the  movement 
of  Joints  1  and  2  is  reversed.  Figures  4.14-4.16  show  how  poorly  the  Tarokh  controller 
responds  to  the  combination  of  different  initial  conditions  and  reversed  trajectories  for 
Joints  1  and  2.  Tarokh’s  controller  is  very  sensitive  to  these  conditions.  Peak  tracking 
errors  for  Joints  1  and  2  are  five  times  larger  for  Trajectory  4  than  for  Trajectory  0. 
Tarokh’s  algorithm  depends  on  position  and  position  error  ternio  only.  Seraji’s  controller 
is  sensitive  to  these  terms  as  well  as  the  desired  velocity  trajectory  components. 
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Figure  4.11.  Effects  of  Tarokh  Discrete  Adaptive  Controller’s  Trajectory  0  Parameters 
on  Trajectoiies  2  and  3 
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Figure  4.12.  Effects  of  Tarokh  Discrete  Adaptive  Controller's  Trajectory  0  Parameters 
on  Tiajectorie&  2  and  3 
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The  movement  of  Trajectory  6  is  the  same  as  Trajectory  1  except  Joint  3  moves  in  the 
opposite  direction.  Desired  velocity  and  acceleration  terms  for  Joint  3  should  have  opposite 
signs.  Figures  4.17-4.19  illustrate  the  Tarokh  controller’s  inability  to  handle  changes  in 
speed,  coupled  with  reversed  trajectory  for  Joint  3.  Figures  4.20-4.22  demonstrate  the 
effects  of  trajectory  speed  and  payload  on  a  Tarokh  controller.  Trajectory  1  causes  the 
same  angular  position  changes  as  Trajectory  0  but  in  1.5  seconds  instead  of  2.0  seconds. 
The  desired  trajectory  components  will  all  be  changed.  Tarokh’s  algorithm  appears  to  be 
unable  to  provide  trajectory  tracking  control  for  conditions  other  than  the  specific  tuning 
trajectory. 

4*5  Summary 

This  chapter  examined  the  decentralized  digital  adaptive  controllers  of  Tarokh  [22] 
and  Seraji  [20].  They  are  proposed  as  computationally  simple,  fast  controllers.  Each 
controller  is  mathematically  developed  as  an  algorithm  that  uses  its  parameters  to  develop  a 
control  torque  as  the  sum  of  auxiliary,  feedforward,  and  feedback  torques.  Both  algorithms 
are  heavily  influenced  by  auxiliary  torque  terms.  In  the  tuning  process,  both  controllers 
tend  to  minimize  the  need  for  one  of  the  other  torques.  Part  of  the  simplicity  of  of  these 
controllers  is  their  independence  from  the  robot’s  dynamics  models,  which  makes  Tarokh’s 
controller  exceedingly  hard  to  tune  and  trajectory  dependent. 

These  controllers  did  not  perform  well  enough  in  this  study  to  be  considered  as  an 
algorithm  suitable  for  gross  motion  control  but,  the  inclusion  of  the  above  mentioned 
auxiliary  torques  into  Model  Based  Adaptive  Input  Controllers,  as  discussed  in  the  next 
chapter,  is  promising. 
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Figure  4.14.  The  Effects  of  Trajectory  0  Tuning  Parameters  on  Trajectories  4  and  5 
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Figure  4.15.  The  Effects  of  Trajectory  0  Tuning  Parameters  on  Trajectories  4  and  5 
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Figure  4.17.  The 
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Figure  4.18.  The  Effects  of  Trajectory  0  Tuning  Parameters  on  Trajectories  1  and  6 
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Figure  4.20.  The  effects  on  Trajectory  0  Tuning  Parameters  on  Speed  -  2  sec  vs  1.5  sec 


Joint  2  Tracking  Errors 


Figure  4.21.  The  effects  on  Trajectory  0  Tuning  Parameters  on  Speed  -  2  sec  vs  1.5  sec 
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Figure  4.22.  The  effects  on  Trajectory  0  Tuning  Parameters  on  Speed  -  2  sec  vs  1.5  sec 
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V.  Model-Based,  Auxiliary  Input  Controllers 


5.1  Model-Based  Auxiliary  Input  Controller  Implementation 

Leahy  used  the  auxiliary  torque  component  of  Seraji’s  controller  and  formed  a  Model- 
Based  Auxiliary  Input  Controller  (MBAIC).  This  controller  couples  the  auxiliary  torque 
with  the  feedforward  and  feedback  torques  developed  by  an  SMBC  controller.  The  MBAIC 
is  expressed  by: 

r  =  Tax  +  T/b  +  Tff  (5.1) 

The  feedforward  torque  term  is  based  on  compensation  of  all  manipulator  system  dynamics. 
The  feedback  gains  are  set  to  the  values  listed  in  Chapter  3.  The  Sera'ii  and  Tarokh 
auxiliary  equations  implemented  are: 

rax  =  f(k)  =  f(k-l)  +  a^{R(k)-hR(k-l)]  +  /3lR(k)-R(k-l)]  (5.2) 
Tax  =  n{k)  =  ??(fc-  l)  +  /3p0e(fc)  +  (/3/-/3p)0e(fc-  1)  (5.3) 

The  role  of  the  auxiliary  input  has  been  reduced  from  eliminating  all  dynamics  disturbances 
to  compensating  for  variations  from  modeled  dynamics  [10].  The  Seraji-based  MBAIC  had 
the  auxiliary  torque  term,  Equation  (5.2),  rewritten  to  allow  adaptation  of  the  matrices 
on  a  joint-by-joint  basis.  A  similar  MBAIC  was  formed  using  Tarokh ’s  auxiliary  torque 
term,  Equation  (5.3).  Tarokh’s  auxiliary  torque  already  adapts  at  the  joint  level.  The 
Tarokh-based  MBAIC  was  retuned  in  the  same  manner  as  described  in  Section  4.3.2.  The 
Seraji-based  MBAIC  was  tested  using  parameters  selected  by  Leahy  and  Seraji  [10,  20]. 
The  auxiliary  torques  are  initialized  to  zero.  The  parameters  are  listed  in  Table  5.1. 


Table  5.1.  MBAIC  Tuning  Parameters  [10] 


Seraji  Based  MBAIC 

Tarokh  Based  MBAIC 

<*=[50,100,100] 

a=[l,l,l] 

/3=[0,0,0] 

=[2000, 2000, 700] 

Wp=[30,40,12] 

Bi={  10,50,5] 

Wv=[20,20,4] 

A,=A2=(-.3,-.3,-.3] 
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5.2  Comparison  of  the  Tarokh  and  Seraji-Based  MBAIC 


Is  there  a  difference  in  the  auxiliary  torques  as  they  are  implemented?  The  auxiliary 
torque  equations  and  the  weighted  errors  can  be  expanded  and  equivalent  terms  compared. 


f(*0 

R  (ft) 


7j(k) 

*.(*) 


f (k  -  1)  +  af[R{k)  +  R {k  -  1)]  +  /3[R(k)  -  R(k  -  1)] 
T  T 

i(k  -  1 )  +  («y  +  /3)R(A)  +  {af  -  (3)R(k  -  1) 

W  pe{k)  +  Wve(k) 

Wpe(k)  +  Wve(k)  ^ 

(W  p+"jr)<k)-^re(k-l) 

rf(k  -  1)  +  pPde(k)  +  (/3j  -  (3 P)6e(k  -  1) 

R2^e(^  —  1)  +  R-3^e(^) 

R2^(&  ""  l)  T  R3^(^) 


Inserting  the  parameters  in  Table  5.1,  the  auxiliary  equations  become: 


f(fc)  =  i(k  —  1)  -f  diag[1003.5, 2010, 4003]e(fc)  +  diag[3.75, 10, 3]e(k  -  1) 

-  diag[1000, 2000, 4000]e(&  -  2) 

rj(k)  =  7}{k  -  1)  +  diag[2180, 2180, 763]e(fc)  -  diag[2277.1, 2233.5, 795.35]e(Jk  -  1) 
+  diag[107.3, 105.3, 37.53]e(fc  -  2) 


Some  of  these  terms  are  harder  than  others  to  select,  i.e.,  R3  is  a  function  of  Tarokh’s 

A  and  a  parameters,  and  Ts  is  not  as  easy  to  manipulate  as  the  other  variables.  Keeping 
/ 

these  two  limitations  in  mind,  the  respective  terms  of  the  weighting  matrices  and  the 
auxiliary  torque  equations  can  be  equated. 


5.3  Model-Based  Auxiliary  Input  Controller  Test  Results 

Both  MBAIC  controllers  were  evaluated  over  all  of  the  test  trajectories.  The  MBAIC 
algorithms  have  a  clear  tracking  advantage  over  their  respective  decentralized  adaptive 


controllers  and  the  SMBC  controller.  Figures  5. 1-5.3  are  representative  of  the  improve¬ 
ments  gained  by  using  the  MBAIC  algorithms  in  place  of  an  SMBC  controller.  The 
plots  in  Appendix  R  contain  the  tracking  errors  produced  by  the  MBAICs  and  compare 
the  MBAIC  to  the  SMBC  used  to  generate  the  feedback  and  feedforward  torques  for  the 
MBAICs.  These  plots  show  that  the  Tarokh-based  MBAIC  controller  is  better  than  or 
comparable  to  the  Seraji-based  model  on  all  test  trajectories.  However,  as  detailed  in 
the  preceding  section,  equivalence  of  these  two  controllers  is  not  easy  to  achieve  in  view 
of  the  origins  of  some  of  the  multipliers  used  in  the  adaptation  of  the  auxiliary  terms. 
The  Seraji-based  MBAIC  may  be  degraded  by  changes  in  the  PUMA.  The  Seraji-based 
MBAIC  was  tuned  only  during  Leahy’s  work  [10],  The  no-load  tracking  errors  are  all  less 
that  3  milliradians.  Trajectory  3  is  the  only  loaded  test  condition  where  errors  exceeded  3 
milliradians. 

5.4  Comparison  of  MBAIC  and  AMBC  Controllers 

Appendix  S  contains  comparisons  of  the  tracking  errors  produced  by  the  MBAIC 
algorithms  and  the  19-parameter  AMBC  controller.  The  AMBC  controller  works  better 
for  all  conditions  except  Trajectory  6,  but  in  several  instances  the  peak  errors  developed 
by  either  trajectory  are  less  that  three  milliradians.  Both  control  methods  provide  a  high 
degree  of  gross  manipulator  control. 

5.5  Summary 

This  chapter  discussed  two  implementations  of  an  MBAIC  controller  that  used  dif¬ 
ferent  auxiliary  torque  equations.  The  similarity  of  the  auxiliary  torque  equations  was 
discussed.  The  tracking  performance  of  these  controllers  was  compared  and  then  com¬ 
pared  to  SMBC  and  AMBC  algorithms. 
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VI.  Conclusions  and  Recommendations 


6.1  Conclusions 

Adaptive  parameter  techniques  have  been  successfully  applied  to  robotic  manipu¬ 
lator  control.  Adaptive  Model-Based  Control  (AMBC)  methods  proved  viable,  as  did 
Model-Based  Auxiliary  Input  Control  (MBAIC)  techniques.  The  decentralized  adaptive 
controllers,  of  which  the  auxiliary  inputs  were  a  component,  were  not  able  to  approximate 
the  degree  of  control  provided  by  AMBC  or  MBAIC  procedures. 

All  three  types  of  controllers  were  experimentally  evaluated  for  the  first  three  links  of 
the  PUMA-560  manipulator.  The  AMBC  method  sought  to  improve  control  by  adapting 
or  calculating  a  set  of  manipulator  parameters.  Our  19-parameter  AMBC  computes  19 
manipulator  parameters  and  leaves  12  at  preselected  values.  When  uninitialized  tests  were 
conducted,  all  31  parameters  were  set  to  zero  and  only  the  19  adapting  parameters  ever 
contributed  to  developing  the  model-based  feedforward  torques.  Initialized  tests  set  all 
31  parameters  to  predetermined  values,  permit  19  parameters  to  change,  and  use  all  31 
parameters  to  calculate  the  feedforward  torque.  Testing  showed  that  adaptation  of  unini¬ 
tialized  parameters  could  not  compare  to  the  results  attained  by  using  properly  initialized 
parameters.  13-  and  16-parameter  models  were  compared  to  the  19-parameter  AMBC. 
The  13-parameter  model  was  heavily  dependent  on  friction  terms  and  did  not  operate  as 
well  as  the  other  options.  The  16-parameter  AMBC  was  already  established  as  an  excellent 
algorithm  for  gross  motion  control.  The  addition  of  three  more  friction  parameters  to  form 
the  19-parameter  AMBC  did  not  improve  upon  the  16-parameter  AMBC  enough  to  justify 
implementation.  The  AMBC  algorithm  is  more  complex  and  computationally  demanding 
than  the  other  two  types  of  controllers  examined. 

Both  decentralized  adaptive  controllers  adapted  six  diagonal  matrices  for  use  in  calcu¬ 
lating  control  torques.  However,  in  practice,  three  of  these  matrices  were  disabled  causing 
Seraji’s  controller  [20]  to  omit,  feedforward  torques  and  the  Tarokh  controller  [22,  23]  t.o 
ignore  part  of  the  feedforward  and  all  of  the  feedback  torques.  Tarokh’s  decentralized 
adaptive  controller  scaled  each  joint’s  adaptation  process  using  diagonal  matrices  of  pro¬ 
portional  and  integral  constants,  a  trajectory  parameter,  and  a  diagonal  weighted  error 
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matrix.  Seraji’s  controller  uses  a  single  diagonal  matrix,  the  trajectory  quantity  used  in 
developing  the  adaptive  matrix’s  torque  component  (i.e.,  C 64  is  a  feedforward  torque  com¬ 
ponent  and  6d  is  used  in  adapting  C  as  in  Equation  (4.9)),  and  diagonal  weighted  error 
and  position  matrix.  The  diagonal  matrix  is  implemented  as  a  scalar.  These  controllers 
have  fewer  characteristics  adapting  on  a  joint  level  than  the  AMBC  controller  and  do  not 
perform  as  well  as  the  AMBC  algorithm. 

The  MBAIC  controllers  use  the  auxiliary  torques  produced  by  the  decentralized  con¬ 
trollers  to  supplement  a  Single  Model-Based  Controller  (SMBC).  This  restores  the  overall 
control  torque  to  its  full  complement  of  feedforward,  feedback,  and  auxiliary  torque  com¬ 
ponents.  These  controllers  easily  outperform  the  decentralized  adaptive  controllers  and 
SMBCs  from  which  they  are  derived.  Our  MBAIC  that  was  based  on  Tarokh’s  auxiliary 
torque  computation  provided  better  tracking  than  the  Seraji  version. 

The  AMBC  controller  worked  better  than  the  MBAIC  algorithms  for  all  test  con¬ 
ditions  except  Trajectory  6.  In  several  instances  the  peak  errors  developed  by  either 
trajectory  are  less  that  three  milliradians.  Both  controllers  provide  a  high  degree  of  gross 
manipulator  control.  The  MBAIC  controllers  are  much  simpler  to  implement  and  compu¬ 
tationally  much  simpler  than  the  AMBC.  If  computational  loading  or  sample  frequencies 
are  a  major  concern,  MBAIC  is  the  method  of  choice.  For  most  applications,  the  AMBC 
controller  should  be  chosen. 

The  aim  of  this  research  was  to  evaluate  control  different  methods  for  implementation 
as  gross  robotic  motion  controllers  to  meet  future  Air  Force  applications.  The  16-parameter 
AMBC  was  verified  as  a  reasonable  approach.  Two  decentralized  adaptive  algorithms  have 
been  rejected  and  the  potential  of  Model-Based  Auxiliary  Input  Control  has  been  identified. 

6.2  Recommendations 

The  decentralized  adaptive  controllers  contain  unnecessary  calculations  in  view  of  the 
number  of  the  mathematically  derived  components  that  were  not  used  in  either  implemen¬ 
tation.  A  sample  period  of  4.5  milliseconds  was  used  for  all  of  the  tests  performed  in  this 
research  effort.  This  software  needs  to  be  rewritten  to  eliminate  avoidable  additions  and 
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multiplications  of  terms  that  will  a 'ways  be  set  to  zero.  These  controllers  may  perform 
admirably  with  an  increased  sampling  rate. 

Additionally,  the  AMBC  and  MBAIC  techniques  used  for  adaptive  robot  control 
should  be  compared  to  other  techniques  proposed  in  control  literature.  This  can  be  accom¬ 
plished  in  the  experimental  environment  available  in  the  Air  Force  Institute  of  Technology 
Robotic  Systems  Laboratory. 
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Appendix  A.  13-Parameter  Uninitialized  Learning  Runs 


This  section  contain  the  results  of  13-parameter  adaptation  runs  for  each  test  tra¬ 
jectory  using  uninitialized  parameters.  Each  trajectory  was  run  nine  times  to  allow  the 
controller  to  adapt. 
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Figure  A.l.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  A. 2. 
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tgure  A. 5.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3 
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Figure  A. 7.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Fayload 
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Figure  A. 8.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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\.9.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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gure  A. 11.  13-Parameter  Uninitialized  Adaptation  Runs  ••  Trajectory  4 

Joint  3  Tracking  Errors 


v/-; 


0.50  1.00  1.50 

Time  (seconds) 


•"N.’-vv 


2.00 


Figure  A. 12.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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gure  A.13.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  5 
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gure  A. 15.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  5 
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Figure  A. 17.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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Figure  A. 19.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  A. 20.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  A. 23.  13-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  6 
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This  section  contain  the  results  of  16-Parameter  adaptation  runs  for  each  test  tra¬ 
jectory  using  uninitialized  parameters.  Each  trajectory  was  run  nine  times  to  allow  the 
controller  to  adapt. 
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igure  B.l.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  B.2.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  B.3.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  B.4.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3 
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Figure  B.6.  16-Paraineter  Uninitialized  Adaptation  Runs  -  Trajeetmy  3 
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Figure  B.7.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  B.8.  1G-Parameter  Uninitialized  Adaptation  Runs  -  Trajectoiy  3  w /  Payload 
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Figure  B.9.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  B.ll.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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Figure  B.12.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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Figure  B.13.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  5 
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Figure  B.14.  16-Parameler  Uninitialized  Adaptation  Runs  ■  Trajectory  5 
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gure  B.15.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  5 
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Figure  B.16.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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B.17.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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Figure  B.18.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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Figure  B.19.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  B.20.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w /  Payload 
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Figure  B.21.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w /  Payload 
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B.23.  16-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  6 
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Appendix  C.  19-Parameter  Uninitialized  Learning  Runs 


This  section  contain  the  results  of  19-parameter  adaptation  runs  for  each  test  tra¬ 
jectory  using  uninitialized  parameters.  Each  trajectory  was  run  nine  times  to  allow  the 
controller  to  adapt. 
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Figure  C.l.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  C.2.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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igure  C.3.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  2 
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Figure  C.4.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3 
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19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3 
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Figure  C.7.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  C.8.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w /  Payload 
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Figure  C.9.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  C.10.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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gure  C.ll.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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Figure  C.12.  19-Paramctcr  Uninitialized  Adaptation  Runs  -  Trajectory  4 
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Figure  C.14.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  5 
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gure  C.17.,  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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Figure  C.18.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1 
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Figure  C.19.  19  Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  C.20.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  C.21.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  C.22.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  G 
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Figure  C.23.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  6 
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Figure  C.24.  19-Parameter  Uninitialized  Adaptation  Runs  -  Trajectory  6 
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Appendix  D.  Comparisons  of  Single  Uninitialized  Runs 


This  section  contains  a  comparison  of  the  first  runs  for  13-,  16-,  and  19-parameter 
AMBC  controllers  for  each  test  trajectory  using  uninitialized  parameters. 
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Figure  D.l.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  2 
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Figure  D.4.  Comparison  of  Single  Uninitia 
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Figure  D.7.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  3  w/  Payload 

Joint  2  Tracking  Errors 


4 


0.035  - 

0.030  - 

0.025  -  i 

0.020  -  | 

0.015  - 

0.010  • 

0.005  - 

0.000  (-• 

-0.005  - 

-0.010  — 

0.00 


0.50  1.00  1.50 

Time  (seconds) 


Figure  D.8.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  3  w/  Payload 
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Figure  D.9.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  3  w/  Payload 
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Figure  D.10.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  4 
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Figure  D.ll.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  4 
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Figure  D.12,  Comparison  of  Single  Uninitia 
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Figure  D.13.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  5 
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Figure  D.14.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  5 
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Figure  D.lb.  Comparison  of  Single  Uninitialized  Huns  -  Trajectory  1 
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Figure  D.17.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  1 
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Figure  D.1S.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  1 
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Figure  D.19.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  1  w/  Payload 
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Figure  D.20.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  l  w/  Payload 
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D.21.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  1  w/  Payload 
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Figure  D.22.  Comparison  of  Single  Uninitialized  Runs  -  Trajectory  6 
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Appendix  E.,  Comparisons  of  Uninitialized  Runs  After  Learning 


This  section  contains  a  comparison  of  the  13- ,  16- ,  and  19-parameter  adaptation  runs 
for  each  test  trajectory  using  uninitialized  parameters.  Each  trajectory  was  run  nine  times 
to  allow  the  controller  to  adapt. 
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Figure  E.2.  Comparison  of  Uninitialized  Runs  After  Learning 
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Figure  E.3.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  2 
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Figure  E.4.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  3 
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Figure  E.5.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  3 
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Figure  E.6.  Comparison  of  Uninitialized  Runs  After  Learning  - 
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Figure  E.7.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  3  w/  Fayload 
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Figure  E.8.  Compcuismi  uf  Uninitialized  Runs  After  Learning  -  Trajectory  3  w /  Payload 
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Figure  E.ll.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  4 
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Figure  E.12.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  4 
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Figure  E.14.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  5 
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Figure  E.15.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  5 
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Figure  E.16.  Comparison  of  Uninitialized  Thins  After  Learning  -  Trajectory  1 
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Figure  E.17.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  1 
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Figure  E.18.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  1 
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Figure  E.19.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  1  w/  Payload 
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Figure  F.9.O.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  1  w /  Payload 
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Figure  R22.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  6 
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Figure  R.24.  Comparison  of  Uninitialized  Runs  After  Learning  -  Trajectory  6 


19  Parameters 

ni§§§§ 

13  Parameters 

16  Parameters 

SMBC 

Appendix  F.  13-Parameter  Initialized  Learning  Runs 


This  section  contain  the  results  of  13-parameter  adaptation  runs  for  each  test  trajec¬ 
tory  using  initialized  parameters.  Each  trajectory  was  run  five  times  to  allow  the  controller 
to  adapt. 
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Figure  F.l«  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  2 
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Figure  F.2.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  2 
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Figure  F.3.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  2 
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Figire  F.4.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3 


First  Run 

Fourth  Run 

IBB 

Second  Run 

bb 

Fifth  Run 

-3 


0.0015 

0.0010 


w  0.0005 


'O  0.0000 
(u 


-0.0005 


M  -0.0010 
W 


-0.0015 


Joint  2  Tracking  Errors 


4  # 


•0.0020 

0.00 


0.50  1.00  1.50 

Time  (seconds) 


Figure  F.5.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3 
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Figure  F.6^13-Parametcr  Initialized  Adaptation  Runs  -  Trajectory  3 
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Figure  F.8.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  F.12.  13-Parameter  Initialized  Adaptation  Puns  - 
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Figure  F.13.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  5 
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Figure  F.15.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  5 
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Figure  F.16.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1 
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Figure  F.18.  13-Parameter  Initialized  Adaptation  Runs  -  1 


First  Run 

Fourth  Run 

Second  Run 

— 

Fifth  Run 

'rajectory  1 


F-10 


Error  (radians)  3  Error  (radians) 


0.004 


Joint  1  Trackin'!  Errors 


Time  (seconds) 

F.19.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 


Joint  2  Tracking  Errors 


Time  (seconds) 


Figure  F.20.  13  Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w /  Payload 
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Figure  F.21.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  F.22.  13-Parameter  Initialized  Adaptation  Runs  -  Trajectory  0 
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Appendix  G.  16-Parameter  Initialized  Learning  Runs 


This  section  contain  the  results  of  16-parameter  adaptation  runs  for  each  test  trajec¬ 
tory  using  initialized  parameters.  Each  trajectory  was  run  five  times  to  allow  the  controller 
to  adapt. 
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Figure  G.3.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  2 
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Figure  G.4.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3 
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Figure  G.7.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3  w/  Payload 
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Figure  G.8.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3  w /  Payload 
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Figure  G.12.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  A 
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Figure  G.13.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  5 
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Figure  G.17.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1 
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Figure  G.19.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  G.20.  16-Parameler  Initialized  Adaptation  Runs  -  Trajectory  1  w /  Payload 
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Figure  G.23.  16-Parameter  Initialized  Adaptation  Runs  -  Trajectory  6 
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Appendix  H.  19-Parameter  Initialized  Learning  Runs 


This  section  contain  the  results  of  19-parameter  adaptation  runs  for  each  test  trajec¬ 
tory  using  initialized  parameters.  Each  trajectory  was  run  five  times  to  allow  the  controller 
to  adapt. 


II-l 


Joint  1  Tracking  Errors 


0.00  0.50  1.00  1.50  2.00 

Time  (seconds) 

Figure  H.l.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  2 


Joint  2  Tracking  Errors 
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Figure  H.5.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3 
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Figure  II.6.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  3 
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figure  H.ll.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  4 
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Figure  H.13.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  5 
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Figure  H.14.  19-Parametcr  Initialized  Adaptation  Runs  -  Trajectory  5 
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■'igure  H.15.  19-Parameter  Initialized  Adaptation  Ituns  -  Trajectory  5 
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Figure  H.16.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1 
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Figure  H.17.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1 
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Figure  H.18.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1 
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H.19.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w/  Payload 
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Figure  H.20.  19-Parameter  Initialized  Adaptation  Runs  -  Trajectory  1  w /  Payload 
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Appendix  I.  Comparisons  of  Single  Initialized  Runs 


This  section  contains  a  comparison  of  the  first  runs  for  13-,  16-,  and  19-parameter 
AMBC  controllers  for  each  test  trajectory  using  initialized  parameters. 


1-1 


Error  (radians)  Error  (radians) 


Joint  1  Tracking  Errors 


Time  (seconds) 

Figure  1.1.  Comparison  of  Single  Initialized  Runs  -  Trajectory  2 


Joint  2  Tracking  Errors 


Figure  T.2.  Comparison  of  Single  Initialized  Huns  -  Trajectory  2 


19  Parameters 

13  Parameters 

16  Parameters 

SMBC 

1-2 


Error  (radians)  Error  (radians) 


Joint  3  Tracking  Errors 


Time  (seconds) 

Figure  1.3. ;  Comparison  of  Single  Initialized  Runs  -  Trajectory  2 


Joint  1  Tracking  Errors 


Figure  1.4.  Comparison  of  Single  Initialized  Runs  -  Trajectory  3 


19  Parameters 

13  Parameters 

16  Parameters 

SMBC 

1-3 


Joint  1  Tracking  Errors 


0.00  0.50  1.00  1.50  2.00 

Time  (seconds) 
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Appendix  J.  Comparisons  of  Initialized  Runs  After  Learning 


This  section  contains  a  comparison  of  the  13-,  16-,  and  19-parameter  adaptation  runs 
for  each  test  trajectory  using  initialized  parameters.  Each  trajectory  was  run  five  times  to 
allow  the  controller  to  adapt. 
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Figure  J.14.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  5 


IB 

19  Parameters 

13  Parameters 

1^5 

16  Parameters 

J-8 


Joint  3  Tracking  Errors 


Time  (seconds) 

Figure  J.15.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  5 
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Figure  J.16.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  1 
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Figure  J.17.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  1 
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Figure  J.18.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  1 
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Figure  J.20.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  1  w/  Payload 
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21.  Comparison  of  Initialized  Runs  After  Learning  -  Trajectory  1  w/  Payload 


Joint  1  Tracking  Errors 


Figure  J.22.  Comparison  of  Initialized  Runs  Alter  Learning  -  Trajectory  6 
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Appendix  K.  Comparisons  of  16-  and  19-Parameter  Learning  Runs 


This  section  contains  a  comparison  of  the  final  learning  runs  for  16-  and  19-parameter 
AMBC  controllers  for  each  test  trajectory  using  uninitialized  and  initialized  parameters. 
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Figure  K.l.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  2 
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Figure  K.2.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  2 
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Figure  K.3.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  2 
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Figure  K.4.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  3 
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Figure  K.6.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  3 
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Figure  K.7.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  3  w/ 
Payload 
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Figure  K.8.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  Trajectory  3  w/ 
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Figure  K.9.  Comparison  of  lo-  and  19-Parameter  Adaptation  Runs  -  Trajectory  3  w/ 
Payload 
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Figure  K.10.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  4 
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Figure  K.ll.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  4 
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Figure  K.12.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  4 
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Figure  K.13.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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Figure  K.14.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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Figure  K.15.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  5 
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Figure  K.16.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  1 
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Figure  K.17.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  1 
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Figure  K.18.  Comparison  of  1G-  and  19-Parameter  Adaptation  Runs  -  Trajectory  1 
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Figure  K.19.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  1  w/ 
Payload 
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Figure  K.20.  Comparison  of  16-  and  19-Paramet.er  Adaptation  Runs  -  Trajectory  1  w / 
Payload 
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Figure  K.21.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  1  w/ 
Payload 
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Figure  K.23.  Comparison  of  16-  and  19-Parameter  Adaptation  Runs  -  Trajectory  6 
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Figure  K.24.  Comparison  of  16  and  19-Paramclcr  Adaptation  Runs  -  Trajectory  6 
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Appendix  L.  Learned  Parameters  for  13-Parameter  Testing 


The  results  of  conducting  learning  runs  for  all  six  test  trajectories  using  13  adaptive 
parameters  are  listed  in  the  following  tables.  These  trajectories  were  run  with  an  assess¬ 
ment  of  all  of  the  manipulator’s  parameters  given  (Initialized)  and  with  all  parameters  set 
to  zero  (w/o  Init).  The  learning  process  was  conducted  over  a  series  of  five  runs  for  the 
initialized  tests  and  nine  runs  for  the  uninitialized  tests. 


Table  L.l.  Learned  13-Parameter  Values  -  Trajectory  1 


Trajectory  1  without  payload  | 

Trajectory  1  with  payload  j 

Number 

Initialized 

w/o  Init 

Nominal 

Initialized 

w/o  Init 

1 

0.30 

0.85 

-4.96 

0.68 

1.37 

-9.75 

2 

-52.09 

-54.86 

-30.73 

-60.56 

-63.99 

-33.77 

3 

-7.53 

-7.30 

-10.12 

-17.12 

-18.59 

-20.13 

4 

-0.01 

-0.01 

0.00 

-0.01 

-0.01 

0.00 

5 

4.50 

-1.73 

-36.53 

4.50 

-2.58 

-38.47 

6 

3.50 

-0.23 

-23.76 

3.50 

3.81 

-19.08 

7 

3.50 

1.09 

-9.25 

3.50 

2.51 

-10.02 

8 

5.95 

7.62 

10.42 

5.95 

7.24 

9.95 

9 

6.82 

10.58 

20.56 

6.82 

8.78 

18.74 

10 

3.91 

4.30 

4.72 

3.91 

1.98 

2.82 

11 

1.00 

0.21 

-5.11 

1.00 

-0.59 

-6.28 

12 

1.00 

-0.34 

-2.43 

1.00 

-0.95 

-5.92 

13 

1.00 

0.37 

-5.90 

1.00 

-1.61 

-11.04 

L-l 


Table  L.2.  Learned  13-Parameter  Values  -  Trajectory  3 


Trajectory  3  without  payload 

Trajectory  3  with  payload 

Number 

Nominal 

Initialized 

w/o  Init 

Nominal 

Initialized 

w/o  Init 

1 

0.30 

-0.16 

-3.28 

0.68 

-0.40 

-7.44 

2 

-52.09 

-52.32 

-15.75 

-60.56 

-61.30 

-19.53 

3 

-7.53 

-8.34 

-10.86 

-17.12 

-19.05 

-24.65 

4 

-0.01 

-0.01 

0.00 

-0.01 

-0.01 

0.00 

5 

4.50 

6.03 

-43.66 

4.50 

6.11 

-50.12 

6 

3.50 

5.91 

-14.67 

3.50 

8.35 

-6.52 

7 

3.50 

2.16 

-8.76 

3.50 

4.37 

-13.15 

8 

5.95 

7.29 

7.81 

5.95 

6.24 

3.30 

9 

6.82 

7.79 

32.64 

6.82 

5.23 

31.98 

10 

3.91 

4.47 

6.84 

3.91 

3.65 

■H 

11 

1.00 

0.69 

-5.05 

1.00 

-0.30 

-4.40 

12 

1.00 

-1.03 

1.75 

1.00 

-2.40 

-4.07 

13 

1.00 

0.00 

-5.61 

1.00 

-2.58 

-7.37 

Table  L.3.  Learned  13-Parameter  Values  -  Trajectories  2  and  4 


Trajectory  2 

Trajectory  4 

Number 

Nominal 

Initialized 

w/o  Init 

Initialized 

w/o  Init 

1 

0.30 

0.26 

-4.66 

0.69 

-4.18 

2 

-52.09 

-51.67 

-8.21 

-52.26 

-24.78 

3 

-7.53 

-8.67 

-9.23 

-7.57 

-7.37 

4 

-0.01 

-0.01 

0.00 

-0.01 

0.00 

5 

4.50 

10.24 

-47.31 

0.51 

-28.24 

6 

3.50 

5.84 

-5.37 

4.07 

-39.84 

7 

3.50 

4.47 

5.33 

2.36 

7.80 

8 

5.95 

6.71 

8.57 

8.11 

7.08 

9 

6.82 

4.27 

47.05 

10.74 

-2.69 

10 

3.91 

5.32 

9.90 

4.75 

7.43 

11 

1.00 

0.76 

-5.45 

0.02 

-0.99 

12 

1.00 

-0.65 

0.54 

-0.15 

-0.44  j 

13 

1.00 

0.19 

-2.92 

0.25 

-2.39 

L-2 


Table  L.4.  Learned  13-Parameter  Values 


Trajectory  5 


Initialized  w/o  Init 


1.25 


-54.46 


-8.17 


-0.01 


0.35 


-1.55 


4.62 


7.71 


4.55 


6.03 


0.28 


-  Trajectories  5  and  6 


6.85 


-32.77 


-11.51 


0.00 


-24.52 


-34.95 


11.56 


8.69 


12.21 


10.87 


-2.44 


Trajectory  6 

Initialized 

w/o  Init 

|  2.29 

4.86 

-46.72 

-36.17 

-10.58 

-6.63 

-0.01 

0.00 

17.40 

-36.41 

21.55 

-43.15 

6.24 

-11.07 

7.96 

7.56 

7.34 

7.20 

5.26 

2.30 

0.55 

-5.07 

-0.48 

-5.30 

0.06 

-7.07 

Appendix  M.  Learned  Parameters  for  16-Parameter  Testing 


The  results  of  conducting  learning  runs  for  all  six  test  trajectories  using  adaptive 
16-parameters  are  listed  in  the  following  tables.  These  trajectories  were  run  with  an 
assessment  of  all  of  the  manipulator’s  parameters  given  (Initialized)  and  with  all  parameters 
set  to  zero  (w/o  Init).  The  learning  process  was  conducted  over  a  series  of  five  runs  for 
the  initialized  tests  and  nine  runs  for  the  uninitialized  tests. 


Table  M.l.  Learned  16-Parameter  Values  -  Trajectory  1 


Trajectory  1  without  payload  | 

Trajectory  1  with  payload 

Number 

Nominal 

Initialized 

w/o  Init 

Nominal 

Initialized 

w/o  Init 

1 

0.30 

1.04 

5.81 

0.68 

0.05 

6.94 

2 

-52.09 

-53.83 

-17.64 

-60.56 

-62.62 

-20.34 

3 

-7.53 

-8.81 

-8.46 

-17.12 

-18.86 

-15.36 

4 

-0.01 

-0.01 

0.00 

-0.01 

-0.01 

0.00 

5 

-0.03 

0.23 

-1.08 

-0.06 

0.34 

-1.40 

6 

-1.57 

-2  86 

-3.34 

-1.46 

-3.40 

-4.33 

7 

2.10 

2.95 

0.65 

2.47 

3.03 

-0.71 

8 

-0.03 

-0.70 

-1.20 

-0.07 

-0.89 

-1.51 

9 

0.67 

2.99 

11.41 

1.51 

2.99 

11.84 

10 

-0.12 

-0.23 

-1.90 

-0.26 

-0.36 

-2.65 

11 

4.50 

4.81 

7.04 

4.50 

5.45 

7.49 

12 

3.50 

4.08 

3.58 

3.50 

4.37 

4.32 

13 

3.50 

3.45 

3.18 

3.50 

3.51 

3.07 

14 

5.95 

7.46 

5.02 

5.95 

6.91 

4.61 

15 

9.51 

17.00 

6.82 

8.50 

18.29 

16 

5.86 

10.14 

3.91 

5.67 

8.30 

M-l 


Table  M.2.  Learned  16-Parameter  Values  -  Trajectory  3 


Trajectory  3  with  payload 


Nominal  Initialized  w/o  Init 


0.68  -0.90  8.65 


-60.56  -61.03  -39.21 


-17.12  -19.03  -21.75 


-0.01  -0.01  0.00 


-0.06  0.83  3.31 


-1.46  -2.50  -3.73 


2.47  2.46 


-0.07  -0.79  1.79 


1.51  1.47  11.79 


-0.26  -0.19  -0.95 


4.50  5.39  7.39 


3.50  4.33  5.16 


3.50  3.70  5.34 


5.95  7.23  7.52 


6.82  6.97  21.92 


3.91  4.66  8.20 


Trajectory  3  without  payload 

Nominal 

Initialized 

w/o  Init 

0.30 

-0.83 

8.70 

-52.09 

-52.29 

-33.41 

-7.53 

-8.85 

-9.02  | 

-0.01 

-0.01 

0.00 

-0.03 

0.87 

4.33 

-1.57 

-2.26 

-2.33 

2.10 

2.13 

6.43 

-0.03 

-0.65 

1.69 

0.67 

0.91 

8.19 

-0.12 

WEH\ 

4.50 

5.26 

5.21 

3.50 

4.15 

3.81 

3.50 

3.55 

3.26 

5.95 

7.84 

10.11 

6.82 

6.88 

20.25 

3.91 

4.48 

11.24 

Table  M.3.  Learned  16-Parameter  Values  -  Trajectories  2  and  4 


Traje  aory  2 

Trajectory  4 

Number 

Nominal 

Initialized 

w/o  Init 

Initialized 

w/o  Init 

i 

0.30 

1.36 

5.46 

0.45 

-3.57 

2 

-52.09 

-51.35 

-29.85 

-50.99 

-16.47 

3 

-7.53 

-9.26 

-17.22 

-7.17 

-3.74 

4 

-0.01 

-0.01 

0.00 

-0.01 

0.00 

5 

-0.03 

-1.03 

-1.03 

0.20 

3.34 

6 

-1.57 

-1.35 

-2.47 

-0.86 

-3.71 

i  7 

2.10 

1.14 

2.94 

2.00 

0.32 

8 

-0.03 

0.11 

1.20 

-0.10 

-1.42 

1  9 

0.67 

0.46 

6.39 

2.11 

10.51 

10 

-0.12 

-0.37 

-0.47 

0.20 

0.69 

11 

4.50 

5.27 

1.93 

5.46 

3.64 

12 

3.50 

3.90 

2.49 

3.93  ! 

-1.95 

j  13 

3.50 

3.56 

0.24 

2.86 

1.82 

14 

5.95 

7.36 

2.27 

8.13 

9.70 

15 

6.82 

6.12 

21.36 

8.86 

-8.96  1 

16 

3.91 

5.40 

4.58 

4.06 

4.63 

M-2 


Table  M.4.  Learned  16-Parameter  Values  -  Trajectories  5  and  6 


Trajectory  5 

Trajectory  6 

Number 

Nominal 

Initialized 

w/o  Init 

Initialized 

w/o  Init 

1 

0.30 

0.40 

-6.43 

-0.29 

-0.56 

2 

-52.09 

-53.11 

-21.81 

-44.49 

-29.35 

3 

-7.53 

-8.14 

-13.48 

-9.92 

6.23 

4 

-0.01 

-0.01 

0.00 

-0.01 

0.00 

5 

-0.03 

0.78 

1.02 

0.24 

2.93 

6 

-1.57 

-1.43 

-1.75 

-0.47 

2.08 

7 

2.10 

1.36 

-6.85 

-2.94 

24.57 

8 

-0.03 

-0.03 

-0.27 

0.02 

1.04 

9 

0.67 

3.65 

15.35 

-1.58 

14.23 

10 

-0.12 

-0.14 

-0.03 

-0.71 

3.05 

11 

4.50 

5.10 

4.18 

7.29 

4.73 

12 

3.50 

3.54 

1.58 

4.55 

-0.97 

13 

3.50 

3.45 

-0.43 

3.68 

5.46 

14 

5.95 

7.00 

3.99 

6.96 

10.59 

15 

6.82 

5.70 

13.34 

9.99 

-3.93 

16 

3.91 

5.73 

5.18 

6.08 

-4.40 
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Appendix  N.  Learned  Parameters  for  19-Parameter  Testing 


The  results  of  conducting  learning  runs  for  all  six  test  trajectories  using  adaptive 
19-parameters  are  listed  in  the  following  tables.  These  trajectories  were  run  with  an 
assessment  of  all  of  the  manipulator’s  parameters  given  (initialized)  and  with  all  parameters 
set  to  zero  (w/o  Init).,  The  learning  process  was  conducted  over  a  series  of  five  runs  for 
the  initialized  tests  and  nine  runs  for  the  uninitialized  tests. 


Table  N.l.  Learned  19-Parameter  Values  -  Trajectory  1 


Trajectory  1  without  payload 

Trajectory  1  with  payload 

Number 

Nominal 

Initialized 

w/o  Init 

Nominal 

Initialized 

w/o  Init 

i 

0.30 

0.87 

6.78 

0.68 

1.31 

6.71 

2 

-52.09 

-53.15 

-17.88 

-60.56 

-62.80 

-20.81 

3 

-7.53 

-9.63 

-8.88 

-17.12 

-18.73 

-15.39 

4 

-0.01 

-0.01 

0.00 

-0.01 

-0.01 

0.00  1 

5 

-0.03 

0.38 

-1.18 

-0.06 

0.11 

-1.49 

6 

-1.57 

-2.77 

■3.36 

-1.46 

-3.45 

-4.41 

7 

2.10 

3.40 

0.82 

2.47 

3.59 

-0.65 

8 

-0.03 

-0.65 

-1.24 

-0.07 

-0.91 

-1.57 

9 

0.67 

3.04 

11.58 

1.51 

3.59 

11.69 

j  10 

-0.12 

-0.16 

-1.80 

-0.26 

-0.25 

-2.62 

11 

4.50 

4.23 

7.00 

4.50 

5.37 

7.57 

12 

3.50 

4.16 

3.61 

3.50 

4.54 

4.44 

13 

3.50 

3.52 

3.33 

3.50 

3.65 

3.20 

i  i4 

5.95 

7.35 

4.74 

5.95 

7.13 

4.54 

15 

6.82 

9.88 

17.13 

6.82 

8.73 

18.04 

16 

3.91 

6.02 

9.45 

3.91 

5.12 

8.21 

17 

1.00 

1.64 

1.22 

1.00 

-0.61 

0.27 

18 

ss 

0.87 

-0.16 

1.00 

0.53 

-0.45 

19 

jhB 

0.91 

-0.33 

1.00 

1.10 

-0.65 

N-l 


Table  N.2.  Learned  19-Parameter  Values  -  Trajectory  3 


Trajectory  3  without  payload 

Trajectory  3  with  payload 

Number 

Nominal 

Initialized 

w/o  Init 

Nominal 

Initialized 

w/o  Init 

1 

0.30 

-1.22 

9.24 

0.68 

-1.24 

8.13 

2 

-52.09 

-51.47 

-33.83 

-60.56 

-60.75 

-38.48 

3 

-7.53 

-7.94 

-10.44 

-17.12 

-19.23 

-18.18 

4 

-0.01 

-0.01 

0.00 

-0.01 

-0.01 

0.00 

5 

-0.03 

1.19 

3.96 

-0.06 

0.70 

4.38 

6 

-1.57 

-2.36 

-2.09 

-1.46 

-3.12 

-4.25  j 

7 

2.10 

2.38 

6.26 

2.47 

2.56 

6.79 

8 

-0.03 

-0.64 

1.66 

-0.07 

-0.98 

1.59 

9 

0.67 

0.64 

8.36 

1.51 

0.86 

10.26 

10 

-0.12 

-0.04 

-0.51 

-0.26 

-0.23 

-1.03 

11 

4.50 

5.10 

5.15 

4.50 

5.70 

6.33 

12 

3.50 

4.30 

3.77 

3.50 

4.70 

5.16 

13 

3.50 

3.56 

3.32 

3.50 

4.11 

4.62 

14 

5.95 

7.61 

10.07 

5.95 

7.11 

9.09 

15 

6.82 

6.42 

20.54 

6.82 

7.01 

21.46 

16 

3.91 

4.65 

10.67 

3.91 

4.88 

9.54 

17 

1.00 

0.44 

-0.63 

1.00 

0.16 

-2.74 

1  18 

1.00 

0.45 

-0.21 

1.00 

0.34 

-1.07 

19 

1.00 

0.80 

0.01 

1.00 

0.50 

-0.83 

N-2 


Table  N.3.  Learned  19-Parameter  Values  -  Trajectories  2  and  4 


Trajectory  2 

Trajectory  4 

Number 

Nominal 

Initialized 

w/o  Init 

Initialized 

w/o  Init 

1 

0.30 

0.52 

6.90 

1.30 

-3.62 

2 

-52.09 

-51.19 

-29.88 

-50.89 

-16.61 

3 

-7.53 

-9.03 

-16.96 

-7.79 

-4.48 

4 

-0.01 

-0.01 

0.00 

-0.01 

0.00 

5 

-0.03 

-0.73 

-1.18 

-0.05 

3.47 

6 

-1.57 

-1.80 

-2.24 

-0.44 

-3.44 

7 

2.10 

1.40 

2.93 

2.82 

0.41 

8 

-0.03 

0.21 

1.13 

-0.06 

-1.56 

9 

0.67 

0.32 

6.64 

1.83 

10.43 

10 

-0.12 

-0.48 

-0.47 

0.20 

0.68 

11 

4.50 

5.05 

2.09 

5.46 

3.80 

12 

3.50 

4.17 

2.72 

4.00 

-1.90 

13 

3.50 

3.08 

0.17 

2.69 

1.56 

14 

5.95 

7.78 

1.79 

8.42 

9.48 

15 

6.82 

5.93 

21.30 

7.60 

-9.34 

16 

3.91 

5.70 

4.89 

4.14 

i  4.43 

17 

1.00 

1.93 

3.96 

1.22 

0.70 

18 

1.00 

0.69 

0.49 

0.75 

-0.19 

19 

1.00 

0.89 

-0.42 

0.96 

-0.38 

N-3 


Table  N.4.  Learned  19-Parameter  Values  -  Trajectories  5  and  6 


Trajectory  5 

Trajectory  6 

Number 

Nominal 

Initialized 

w/o  Init 

Initialized 

w/o  Init 

1 

0.30 

0.81 

-6.71 

0.88 

-0.52 

2 

-52.09 

-53.37 

-21.84 

-42.77 

-29.18 

3 

-7.53 

-8.08 

-13.66 

-0.50 

6.65 

4 

-0.01 

-0.01 

0.00 

-0,01 

0.00 

5 

-0.03 

0.73 

1.09 

-5.17 

2.94 

6 

-1.57 

-1.44 

-1.75 

-0.41 

1.98 

7 

2.10 

0.18 

-6.83 

-2.47 

24.47 

8 

-0.03 

0.03 

-0.21 

0.03 

1.10 

9 

0.67 

3.85 

15.60 

7.18 

14.06 

10 

-0.12 

-0.11 

-0.08 

-0.69 

3.08 

11 

4.50 

4.85 

4.56 

3.43 

KB 

12 

3.50 

3.62 

1.64 

4.75 

-0.92 

13 

3.50 

3.44 

-0.33 

9.66 

5.53 

14 

5.95 

6.97 

3.50 

6.01 

10.61 

15 

6.82 

4.91 

12.94 

1.73 

-3.96 

16 

3.91 

5.50 

5.00 

6.12 

-4.24 

17 

1.00 

2.22 

4.50 

0.73 

-0.78 

18 

1.00 

0.59 

0.01 

0.70 

-0.20 

19 

1.00 

0.68 

-0.08 

0.01 

-0.70 

Appendix  0.  Comparison  of  Learned  Parameters 


The  results  of  conducting  learning  runs  for  all  six  test  trajectories  using  13,  16,  and 
19  adaptive  parameters  are  listed  in  the  following  tables.  These  trajectories  were  run 
with  an  assessment  of  all  of  the  manipulator’s  parameters  given  (Initialized)  and  with  all 
parameters  set  to  zero  (w/o  Init).  The  learning  process  was  conducted  over  a  seiies  of 
five  runs  ior  the  initialized  tests  and  nine  runs  for  the  uninitialized  tests.  The  last  nine 
parameters  in  the  13-parameter  adaptation  are  placed  so  that  they  coincide  with  the  same 
friction  parameters  in  the  16-  and  19-parameter  configurations. 


Table  0.1.  Learned  Parameter  Values  -  Trajectory  1 


Trajectory  1  Initialized 

Trajectory  1  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

i 

0.30 

0.85 

1.04 

0.87 

-4.96 

5.81 

6.78 

2 

-52.09 

-54.86 

-53.83 

-53.15 

-30.73 

-17.64 

-17.88 

3 

-7.53 

-7.30 

-8.81 

-9.63 

-10.12 

-8.46 

-8.88 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.23 

0.38 

-1.08 

-1.18 

6 

-1.57 

-2.86 

-2.77 

-3.34 

-3.36 

7 

2.10 

2.95 

3.40 

0.65 

0.82 

8 

-0.03 

-0.70 

-0.65 

-1.20 

-1.24 

9 

0.67 

2.99 

3.04 

11.41 

11.58 

10 

-0.12 

-0.23 

-0.16 

-1.90 

-1.80 

11 

4.50 

-1.73 

4.81 

4.23 

-36.53 

7.04 

7.00 

12 

3.50 

-0.23 

4.08 

4.16 

-23.76 

3.58 

3.61  i 

13 

3.50 

1.09 

3.45 

3.52 

-9.25 

3.18 

3.33 

14 

5.95 

7.62 

7.46 

7.35 

10.42 

5.02 

4.74 

15 

6.82 

10.58 

9.51 

9.88 

20.56 

17.00 

17.13 

16 

3.91 

4.30 

5.86 

6.02 

4.72 

10.14 

9.45 

17 

1.00 

0.21 

1.64 

-5.11 

1.22 

18 

1.00 

-0.34 

0.87 

-2.43 

-0.16 

19 

1.00 

0.37 

0.91 

-5.90 

-0.33 

0-1 


Table  0.2.  Learned  Parameter  Values  -  Trajectory  1  with  payload 


L 

Trajectory  1  Initialized 

Trajectory  1  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.68 

1.37 

0.05 

1.31 

-9.75 

6.94 

6.71 

2 

-60.56 

-63.99 

-62.62 

-62.80 

-33.77 

-20.34 

-20.81 

3 

-17.12 

-18.59 

-18.86 

-18.73 

-20.13 

-15.36 

-15.39 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.06 

0.34 

0.11 

-1.40 

-1.49 

6 

-1.46 

-3.40 

-3.45 

-4.33 

-4.41 

7 

2.47 

3.03 

3.59 

-0.71 

-0.65  | 

8 

-0.07 

-0.89 

-0.91 

-1.51 

-1.57 

9 

1.51 

2.99 

3.59 

11.84 

11.69 

10 

-0.26 

-0.36 

-0.25 

-2.65 

-2.62 

11 

4.50 

-2.58 

5.45 

5.37 

-38.47 

7.49 

7.57 

12 

3.50 

3.81 

4.37 

4.54 

-19.08 

4.32 

4.44 

13 

3.50 

2.51 

3.51 

3.65 

-10.02 

3.07 

3.20 

14 

5.95 

7.24  | 

6.91 

7.13 

9.95 

4.61 

4.54 

15 

6.82 

8.78 

8.50 

8.73 

18.74 

18.29 

18.04 

16 

3.91 

1  -  . 

1.98 

5.67 

5.12 

2.82 

8.30 

8.21 

17 

1.00 

-0.59 

-0.61 

-6.28 

0.27 

18 

1.00 

-0.95 

0.53 

-5.92 

-0.45 

19 

1.00 

-1.61 

1.10 

-11.04 

-0.65 

0-2 


Table  0.3.  Learned  Parameter  Values  -  Trajectory  2 


Trajectory  2  Initialized 

Trajectory  2  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.30 

0.26 

1.36 

0.52 

-4.66 

5.46 

6.90 

2 

-52.09 

-51.67 

-51.35 

-51.19 

-8.21 

-29.85 

-29.88 

3 

-7.53 

-8.67 

-9.26 

-9.03 

-9.23 

-17.22 

-16.96 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

-1.03 

-0.73 

-1.03 

-1.18 

6 

-1.57 

-1.35 

-1.80 

-2.47 

-2.24 

7 

2.10 

1.14 

1.40 

2.94 

2.93 

8 

-0.03 

0.11 

0.21 

1.20 

1.13 

9 

0.67 

0.46 

0.32 

6.39 

6.64 

10 

-0.12 

-0.37 

-0.48 

-0.47 

-0.47 

11 

4.50 

10.24 

5.27 

5.05 

-47.31 

1.93 

2.09 

12 

3.50 

5.84 

3.90 

4.17 

-5.37 

2.49 

2.72 

13 

j  3.50 

4.47 

3.56 

3.08 

5.33 

0.24 

0.17 

14 

5.95 

6.71 

7.36 

7.78 

8.57 

2.27 

1.79 

15 

6.82 

4.27 

6.12 

5.93 

47.05 

21.36 

21.30 

16 

3.91 

5.32 

5.40 

5.70 

9.90 

4.58 

4.89 

17 

1.00 

0.76 

1.93 

-5.45 

3.96 

18 

1.00 

-0.65 

0.69 

0.54 

0.49 

19 

1.00 

0.19 

0.89 

-2.92 

-0.42 

0-3 


Table  0.4.  Learned  Parameter  Values  -  Trajectory  3 


Trajectory  3  Initialized 

Trajectory  3  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1  1 

0.30 

-0.16 

-0.83 

-1.22 

-3.28 

8.70 

9.24 

2 

-52.09 

-52.32 

-52.29 

-51.47 

-15.75 

-33.41 

-33.83 

3 

-7.53 

-8.34 

-8.85 

-7.94 

-10.86 

-9.02 

-10.44 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.87 

1.19 

4.33 

3.96 

6 

-1.57 

-2.26 

-2.36 

-2.33 

-2.09 

7 

2.10 

2.13 

2.38 

6.43 

6.26 

8 

-0.03 

-0.65 

-0.64 

1.69 

1.66 

9 

0.67 

0.91 

0.64 

8.19 

8.36 

10 

-0.12 

-0.05 

-0.04 

-0.73 

-0.51 

11 

4.50 

6.03 

5.26 

5.10 

-43.66 

5.21 

5.15 

12 

3.50 

5.91 

4.15 

4.30 

-14.67 

3.81 

3.77 

13 

3.50 

2.16 

3.55 

3.56 

-8.76 

3.26 

3.32 

14 

5.95 

7.29 

7.84 

7.61 

7.81 

10.11 

10.07 

15 

6.82 

7.79 

6.88 

6.42 

32.64 

20.25 

20.54 

1G 

3.91 

4.47 

4.48 

4.65 

6.84 

11.24 

10.67 

17 

1.00 

0.69 

0.44 

-5.05 

-0.63 

18 

1.00 

-1.03 

0.45 

1.75 

-0.21 

19 

1.00 

0.00 

0.80 

-5.61 

0.01 

0-4 


Table  0.5.  Learned  Parameter  Values  -  Trajectory  3  with  payload 


Trajectory  3  Initialized 

Trajectory  3  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.68 

-0.40 

-0.90 

-1.24 

-7.44 

8.65 

8.13 

2 

-60.56 

-61.30 

-61.03 

-60.75 

-19.53 

-39.21 

-38.48 

3 

-17.12 

-19.05 

-19.03 

-19.23 

-24.65 

-21.75 

-18.18 

j  4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.06 

0.83 

0.70 

3.31 

4.38 

6 

-1.46 

-2.50 

-3.12 

-3.73 

-4.25 

7 

2.47 

2.46 

2.56 

7.08 

6.79 

j  8 

-0.07 

-0.79 

-0.98 

1.79 

1.59 

9 

1.51 

1.47 

0.86 

11.79 

10.26 

10 

-0.26 

-0.19 

-0.23 

-0.95 

-1.03 

11 

4.50 

6.11 

5.39 

5.70 

-50.12 

7.39 

6.33 

12 

3.50 

8.35 

4.33 

4.70 

-6.52 

5.16 

5.16 

13 

3.50 

4.37 

3.70 

4.11 

-13.15 

5.34 

4.62 

14 

5.95 

6.24 

7.23 

7.11 

3.30 

7.52 

9.09 

15 

6.82 

5.23 

6.97 

7.01 

31.98 

21.92 

21.46 

16 

3.91 

3.65 

4.66 

4.88 

1.98 

8.20 

9.54 

17 

1.00 

-0.30 

0.16 

-4.40 

-2.74 

18 

1.00 

-2.40 

0.34 

-4.07 

-1.07 

19 

1.00 

-2.58 

0.50 

-7.37 

-0.83 

0-5 


Table  0.6.  Learned  Parameter  Values  -  Trajectory  4 


Trajectory  4  Initialized 

Trajectory  4  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.30 

0.69 

0.45 

1.30 

-4.18 

-3.57 

-3.62 

2 

-52.09 

-52.26 

-50.99 

-50.89 

-24.78 

-16.47 

-16.61 

3 

-7.53 

-7.57 

-7.17 

-7.79 

-7.37 

-3.74 

-4.48 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.20 

-0.05 

3.34 

3.47 

6 

-1.57 

-0.86 

-0.44 

-3.71 

-3.44 

7 

2.10 

2.00 

2.82 

0.32 

0.41 

8 

-0.03 

-0.10 

-0.06 

-1.42 

-1.56 

9 

0.67 

2.11 

1.83 

10.51 

10.43 

10 

-0.12 

0.20 

0.20 

0.69 

0.68 

11 

4.50 

0.51 

5.46 

5.46 

-28.24 

3.64 

3.80 

12 

3.50 

4.07 

3.93 

4.00 

-39.84 

-1.95 

-1.90 

13 

3.50 

2.36 

2.86 

2.69 

7.80 

1.82 

1.56 

14 

5.95 

8.11 

8.13 

8.42 

7.08 

9.70 

9.48 

|  15 

6.82 

10.74 

8.86 

7.60 

-2.69 

-8.96 

-9.34 

16 

3.91 

4.75 

4.06 

4.14 

7.43 

4.63 

4.43 

17 

1  00 

0.02 

-0.99 

0.70 

18 

-0.15 

0.75 

-0.44 

-0.19 

19 

1.00 

0.25 

0.96 

-2.39 

-0.38 

0-6 


Table  0.7.  Learned  Parameter  Values  -  Trajectory  5 


Trajectory  5  Initialized 

Trajectory  5  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.30 

1.25 

0.40 

0.81 

6.85 

-6.43 

-6.71 

2 

-52.09 

-54.46 

-53.11 

-53.37 

-32.77 

-21.18 

-21.84 

3 

-7.53 

-8.17 

-8.14 

-8.08 

-11.51 

-13.48 

-13.66 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.78 

0.73 

1.02 

1.09 

6 

-1.57 

-1.43 

-1.44 

-1.75 

-1.75 

7 

2.10 

1.36 

0.18 

-6.85 

-6.83 

8 

-0.03 

-0.03 

0.03 

-0.27 

-0.21 

9 

0.67 

3.65 

3.85 

15.35 

15.60 

10 

-0.12 

-0.14 

-0.11 

-0.03 

-0.08 

11 

4.50 

0.35 

5.10 

4.85 

-24.52 

4.18 

4.56 

12 

3.50 

-1.55 

3.54 

3.62 

-34.95 

1.58 

1.64 

13 

3.50 

4.62 

3.45 

3.44 

11.56 

-0.43 

-0.33 

14 

5.95 

7.71 

7.00 

6.97 

8.69 

3.99 

3.50 

15 

6.82 

4.55 

5.70 

4.91 

12.21 

13.34 

12.94 

16 

3.91 

6.03 

5.73 

5.50 

10.87 

5.18 

5.00 

17 

1.00 

0.28 

i 

i 

2.22 

-2.44 

4.50 

18 

1.00 

0.07 

fiWifeU 

\MEM 

19 

1.00 

0.00 

0.68 

-2.82 

— 

0-7 


Table  0.8.  Learned  Parameter  Values  -  Trajectory  6 


Trajectory  6  Initialized 

Trajectory  6  Uninitialized 

Number 

Nominal 

13  par 

16  par 

19  par 

13  par 

16  par 

19  par 

1 

0.30 

2.29 

-0.29 

0.88 

4.86 

-0.56 

-0.52 

2 

-52.09 

-46.72 

-44.49 

-42.77 

-36.17 

-29.35 

-29.18 

3 

-7.53 

-10.58 

-9.92 

-0.50 

-6.63 

6.23 

6.65 

4 

-0.01 

-0.01 

-0.01 

-0.01 

0.00 

0.00 

0.00 

5 

-0.03 

0.24 

-5.17 

2.93 

2.94 

6 

-1.57 

-0.47 

-0.41 

2.08 

1.98 

7 

2.10 

-2.94 

-2.47 

24.57 

24.47 

8 

-0.03 

0.02 

0.03 

1.04 

1.10 

9 

0.67 

-1.58 

7.18 

14.23 

14.06 

10 

-0.12 

-0.71 

-0.69 

3.05 

3.08 

11 

4.50 

17.40 

7.29 

3.43 

-36.41 

4.73 

■^1 

12 

3.50 

21.55 

4.55 

4.75 

-43.15 

-0.97 

-0.92 

13 

3.50 

6.24 

3.68 

9.66 

-11.07 

5.46 

5.53 

14 

5.95 

7.96 

6.96 

6.01 

7.56 

10.59 

10.61 

15 

6.82 

7.34 

9.99 

1.73 

7.20 

-3.93 

-3.96 

16 

3.91 

5.26 

6.08 

6.12 

2.30 

-4.40 

-4.24 

17 

1.00 

0.55 

0.73 

-5.07 

-0.78 

18 

1.00 

-0.48 

0.70 

-5.30 

-0.20 

19 

1.00 

0.06 

0.01 

-7.07 

-0.70 

0-8 


Appendix  P.  Tarokh  Algorithm  Runs 


This  section  contains  plots  comparing  the  errors  produced  on  Trajectory  0  versus 
other  test  trajectories  using  the  set  of  design  parameters  chosen  for  Trajectory  0. 
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Figure  P.4.  The 
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Figure  P.7.  Payload  Effects  on  Trajectory  0  Tuning  Parameters 
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Figure  P.ll.  The  Effects  of  Trajectory  0  Tuning  Parameters  on  Trajectories  1  and  6 
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Figure  P.13.  The  effects  on  Trajectory  0  Tuning  Parameters  on  Speed  -  2  sec  vs  1.5  sec 
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Figure  P.14.  The  effects  on  Trajectory  0  Tuning  Parameters  on  Speed  -  2  sec  vs  1.5  sec 
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Appendix  Q.  Comparisons  of  Tarokh  and  Seraji  Control  Algorithms 


This  appendix  contains  comparison  plots  of  Tarokh’s  and  Seraji’s  adaptive  control 
algorithms  for  all  test  trajectories. 
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Figure  Q.5.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  0  w /  Payload 
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Figure  Q.6.  Comparison  cl  Tarokh  and  Seraji  Controllers  -  Trajectory  0  w/  Payload 
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gure  Q.7.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  2 
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Figure  Q.8.  Comparison  of  Tarokh  and  Seraji  Controlleis  -  Trajectory  2 


Tarokh  Controller 


Seraji  Controller 


Q-5 


Joint  3  Tracking  Errors 


Figure  Q.9.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  2 
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Figure  Q.10.  Comparison  of  Tarokh  and  Seiaji  Conti  ulleib  -  Tiajectoiy  3 
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Figure  Q.13.  Comparison  of  Tarckh  and  Seraji  Controllers  -  Trajectory  3  w/  Payload 
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Figure  Q.14.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  3  w/  Payload 
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•.15.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  3  w/  Payload 
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Figure  Q.16.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  4 


Tarokh  Controller 


Seraji  Controller 


Q-9 


Joint  2  Tracking  Errors 


Figure  Q.17.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  4 
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Figure  Q.18.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  4 
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Figure  Q.19.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  5 
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Figure  Q.20.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  5 
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Figure  Q.21.;  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  5 
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Figure  Q.22.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  1 
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Figure  Q.23.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  1 
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Figure  Q.24.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  1 
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Figure  Q.25.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  1  w/  Payload 
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Figure  Q.26.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  1  w/  Payload 
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).28.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  6 
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Figure  Q.29.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  6 
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Figure  Q.30.  Comparison  of  Tarokh  and  Seraji  Controllers  -  Trajectory  6 
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Appendix  R.  Comparisons  of  MB AlC  and  SMBC  Controllers 


This  appendix  contains  comparison  plots  of  the  tracking  errors  generated  by  the 
Model-Based  Auxiliary  Input  Controllers  (MBAIC)  constructed  from  the  auxiliary  torque 
components  of  Tarokh’s  and  Seraji’s  adaptive  control  algorithms.  These  plots  illustrate 
the  improvements  the  MBAIC  controllers  have  over  the  SMBC  controller  composed  of  only 
PD  feedback  torques  and  a  model-based  feedforward  torque. 
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Figure  R.l.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0 
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Figure  R.2.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0 
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Figure  R.3.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0 
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Figure  R.4.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0  w/  payload 
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.5.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0  w/  Payload 
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Figure  R.6.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  0  w/  Payload 
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Figure  R.8.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  2 
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Figure  R.9.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  2 


Joint  1  Tracking  Errors 


Figure  R.10.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  3 
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Figure  R.15.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  3  w/  Payload 
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Figure  R.16.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  4 
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Figure  R.17.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  4 
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Figure  R.18.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  4 
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Figure  R.21.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  5 
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Figure  R.22.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  1 
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Figure  R.24.  Comparison  of  MBAIC  and  SMBC  Controllers  Trajectory  1 
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Figure  R.25.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  1  w/  Payload 
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27.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  1  w/  Payload 
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Figure  11.28.  Comparison  of  MBAIC  and  SMBC  Controllers  Trajectory  6 
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Figure  R..30.  Comparison  of  MBAIC  and  SMBC  Controllers  -  Trajectory  6 
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Appendix  S.  Comparison  of  MBAIC  and  AMBC  Controllers 

This  appendix  contains  plots  which  compare  the  results  achieved  in  19-parameter 
AMBC  algorithm  testing  with  the  trajectory  errors  produced  in  MBAIC  testing. 
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Figure  S.2.  Comparison  of  MBAIC  and  AMBC  Controllers  Trajectory  2 
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Figure  S.3.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  2 
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Figure  S.d.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  3 
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Figure  S.6.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  3 
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Figure  S.7.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  3  w/  Payload 
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Figure  S.8.  Comparison  of  MBAIC  and  AMBC  Controllers  Trajectory  3  w /  Payload 
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Figure  S.9.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  3  w /  Payload 
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Figure  S.10.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  4 
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Figure  S.ll.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  4 
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Figure  S.12.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  4 
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Figure  S.13.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  5 
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Figure  S.14.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  5 
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Figure  S.15.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  5 
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Figure  S.16.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  1 
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Figure  S.17.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  1 
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Figure  S.18.  Comparison  of  MBAIC  and  AMBC  Controllers  Trajectory  1 
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Figure  S.19.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  1  w/  Payload 


Joint  2  Tracking  Errors 


0.00  0.50  1.00  1.50 

Time  (seconds) 

Figure  S.20.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  1  w/  Payload 
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21.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  1  w/  Payload 
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Figure  S.22.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  6 
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Figure  S.23.  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  6 
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Figure  S.24,  Comparison  of  MBAIC  and  AMBC  Controllers  -  Trajectory  6 
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